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PREFACE TO THE THIRD EDITION 



Many Earth science departments today recruit signifi- 
cant numbers of environmental geoscience students 
alongside mainstream geology or Earth science stud- 
ents. This student-led change in academic emphasis 
suggests that a successor to Chemical Fundamentals of 
Geology will meet the needs of students better if it cov- 
ers the chemical foundations of environmental Earth 
science as well as geology The shift in the emphasis of 
this new edition is reflected in its expanded title. 

This change in emphasis has mostly been accomm- 
odated by expanding the existing chapters of the book. 
Reviewers have, however, highlighted the thin cover- 
age of isotopes in earlier editions and - since isotopes 
contribute so much to our understanding of processes 
within the Earth and on its surface - 1 have devoted a 
new chapter (Ch. 10) to simple concepts of isotope geo- 
chemistry. Although to many students this chapter 
may seem more challenging - conceptually and math- 
ematically - than earlier chapters in the book, I hope it 
will nonetheless help them to appreciate the unique 
insights that isotopes contribute to our science. 



This book was originally conceived to support Earth 
science students having no background at all in chem- 
istry. Despite this intention, teachers using the book 
continue to mention students who "feel that it is too 
high-level for their needs". With the requirements of 
those students in mind, I have tried to lower still fur- 
ther the entry hurdles that the book places in their 
path, for example, by introducing boxes like Box 1.1 
'What is energy?' 

As in previous editions (despite criticism in one pub- 
lished review), where variables and units appear on 
graph axes they are separated by a forward slash (i.e. 
variable / unit) according to Royal Society convention. 1 
The logic of this implied division (variable 4- unit) is 
that, when distance x is divided by a metre, say, or time 
t is divided by a second, the result is a dimensionless 
pure number - numbers (not quantities) are after all 
what we actually plot on graph axes. 



1 See also www.animations.physics.unsw.edu.au/jw/graphs. 
htm#Units. 



PREFACE TO THE SECOND EDITION 



Since its first publication in 1989, this book has been 
used by a much wider range of students than I orig- 
inally envisaged. My initial objective was simply to 
equip elementary Earth science undergraduates with 
the basic chemical understanding that their geological 
degree programme would require, but in the course 
of writing the first edition I fell into the trap of trying 
to entice them a little further so that they would 
embrace some of the underlying scientific principles. 
Accordingly, though somewhat to my surprise, the 
book has been adopted at the advanced undergrad- 
uate level and even as remedial reading for postgrad- 
uate studies as well as for its original purpose. It has 
also been adopted in some chemistry departments as 
a geochemistry primer. 

In preparing this edition I have tried to increase the 
book's usefulness in two directions. Firstly, I have 
sought to make it more reader-friendly to those, such 
as mature students, who may remember very little 
chemistry from their school days or whose experiences 
of school chemistry were discouraging. Secondly, with 



more advanced students in mind, I have introduced a 
little more chemical rigour and more pointers to geo- 
logical and geochemical applications, and have, I hope, 
also reflected some recent advances in cosmological, 
geological and environmental understanding. The 
change in format has allowed me to make more effec- 
tive use of the boxes and to enhance some figures. 

I am grateful to many people who have fed back to 
me their reactions to the first edition, particularly Paul 
Browning, Hilary Downes, Mike Henderson, Bob 
Major, Steven Richardson and Andy Saunders. Not 
every suggestion, of course, has seemed to me consist- 
ent with the original spirit of the book, but all com- 
ments have been valuable in giving me a reader's or a 
teacher's eye view. The encouragement I have received 
from colleagues such as Derek Blundell and Euan 
Nisbet has also been greatly appreciated. My family 
has been as understanding (or as resigned) as ever! 

I thank Lynne Blything for her kind help with figure 
revision, and Ruth Cripwell and Ian Francis for their 
editorial advice and encouragement. 



PREFACE TO THE FIRST EDITION 



Chemical principles are fundamental to a large area of 
geological science, and the student who reads for a 
degree in geology without knowing any basic chemis- 
try is handicapped from the start. This book has been 
written with such students in mind, but I hope it will 
also provide a helpful refresher course and useful 
background reading in 'geo-relevant' chemistry for all 
earth science students who believe in understanding 
rather than merely memorizing. 

The book has been conceived in three broad sections. 
The first. Chapters 1-4, deals with the basic physical 
chemistry of geological processes, emphasizing how 
consideration of energy can aid our understanding. 
The second section. Chapters 5-8, introduces the wave- 
mechanical view of the atom and describes from that 
perspective the various types of chemical bonding that 
give minerals their distinctive properties. The final sec- 
tion, Chapters 9 and 10, surveys the geologically rel- 
evant elements, concluding with a chapter on why 
some are more abundant than others, in the universe 
as a whole and in the Earth in particular. The emphasis 
throughout is on geological and environmental rel- 
evance; laboratory reactions are hardly mentioned. 



The book is designed to be accessible and stimulat- 
ing, even to students who remember none of their 
school chemistry; for these readers, introductory back- 
ground material and a glossary are provided in the 
Appendices. The more advanced material, which I 
hope will sustain the interest of the chemically literate 
reader, has been segregated into boxes that can be 
ignored on a first reading. The text is thoroughly cross- 
referenced to help the browser or the searcher to find 
what each wants. The book might also assist teachers 
of more advanced courses in geochemistry, mineral- 
ogy and petrology by relieving them of the need to 
cover elementary material in class. 

Many colleagues have given me advice and encour- 
agement during the writing of the book, and I am par- 
ticularly grateful to the following colleagues who have 
read and commented on individual chapters and sug- 
gested innumerable improvements: David Alderton, 
Peter Barnard, Keith Cox, Giles Droop, Paul Henderson, 
Steve Killops, Robert Hutchison, Philip Lee and Eric 
Whittaker. Professor W.D. Carlson, Dr T.K. Halstead 
and Dr J.B. Wright read the whole manuscript and pro- 
vided a wealth of helpful comment. The errors and 
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PREFACE TO THE FIRST EDITION 



eccentricities that remain are of course my responsibil- 
ity alone. The book has grown from a lecture course I 
gave at the former Geology Department of Chelsea 
College, and I warmly acknowledge the opportunities 
which that department provided. 

I would like to thank Joan Hirons, Sue Clay and 
Jennifer Callard for typing some of the chapters, and 
Neil Holloway and Christine Flood who prepared a 



number of the figures. I owe a great deal to Roger Jones 
of Unwin Hyman, whose faith in the project sustained 
me in the face of initial difficulties. 

The biggest debt I owe is to Mary, Joanna and Tim, 
who have tolerated with remarkably good humour my 
neglect of family activities, the perpetual retreat up to 
the study, and the curtailment of holidays. I promise 
not to write another book for a long time. 
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1 ENERGY IN GEOCHEMICAL 
PROCESSES 



Introduction 



The purpose of this book is to introduce the average 
Earth science student to chemical principles that are 
fundamental to the sciences of geology and environ- 
mental geoscience. There can be no more fundamental 
place to begin than with the topic of energy (Box 1.1), 
which lies at the heart of both geology and chemistry. 
Energy plays a role in every geological process, from 
the atom-by-atom growth of a mineral crystal to the 
elevation and subsequent erosion of entire mountain 
chains. Consideration of energy provides an incisive 
intellectual tool for analysing the workings of the com- 
plex geological world, making it possible to extract 
from this complexity a few simple underlying princi- 
ples upon which an orderly understanding of Earth 
processes can be based. 

Many natural processes involve a flow of energy. 
The spontaneous melting of an ice crystal, for example. 



requires energy to be transferred by means of heat into 
the crystal from the 'surroundings' (the air or water 
surrounding the crystal). The crystal experiences an 
increase in its internal energy, which transforms it 
into liquid water. The process can be symbolized by 
writing down a formal reaction: 

H 2 0^H 2 0 (1.1) 

ice water 

in which molecules of water (H,0) are represented as 
migrating from the solid state (left-hand side) into the 
liquid state (right-hand side). Even at 0°C, ice and 
water both possess internal energy associated with the 
individual motions of their constituent atoms and mol- 
ecules. This energy content, which we loosely visualize 
as heat 'stored' in each of the substances, is more cor- 
rectly called the enthalpy 1 (symbol H). Because the 

1 Words in bold type indicate terms that are defined in the 
Glossary. 
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CHAPTER 1 



f 

Box 1.1 What is energy? 



The concept of energy is fundamental to all branches of 
science, yet to many people the meaning of the term 
remains elusive. In everyday usage it has many shades of 
meaning, from the personal to the physical to the mystical. 
Its scientific meaning, on the other hand, is very precise. 

To understand what a scientist means by energy, 
the best place to begin is with a related - but more 
tangible - scientific concept that we call work. Work is 
defined most simply as motion against an opposing force 
(Atkins, 2010, p. 23). Work is done, for example, when a 
heavy object is lifted a certain distance above the ground 
against the force of gravity (Figure 1.1.1). The amount of 
work this involves will clearly depend upon how heavy 
the object is, the vertical distance through which its 
centre of gravity is lifted (Figure 1.1.1b), and the strength 
of the gravitational field acting on the object. The work 
done in this operation can be calculated using a simple 
formula: 

\Notk = mxhxg (111) 

J kg; m ms 2 

where m represents the mass of the object (in kg), h is 
the distance through which its centre of gravity is raised 
(in m - see footnote) 2 , and g, known as the acceleration 
due to gravity (metres per second per second = ms~ 2 ), is a 
measure of the strength of the gravitational field where 
the experiment is being carried out; at the Earth’s sur- 
face, the value of g is 9.81 ms 2 . The scientific unit that 
we use to measure work is called the joule (J), which as 



2 m in italics represents mass (a variable in this equation); m in 
regular type is the abbreviation for metres (units). 

V 



Equation 1.1.1 shows is equivalent to kgxmxms^ 2 = 
kgm 2 s~ 2 (see Table A2, Appendix A). Alternative forms of 
work, such as cycling along a road against a strong oppos- 
ing wind, or passing an electric current through a resistor, 
can be quantified using comparably simple equations, but 
whichever equation we use, work is always expressed in 
joules. 

The weight suspended in its elevated position 
(Figure 1.1.1b) can itself do work. When connected to suit- 
able apparatus and allowed to fall, it could drive a pile into 
the ground (this is how a pile-driver works), hammer a nail 
into a piece of wood, or generate electricity (by driving a 
dynamo) to illuminate a light bulb. The work ideally recov- 
erable from the elevated weight in these circumstances is 
given by Equation 1.1.1. If we were to raise the object 
twice as far above the ground (Figure 1.1.1c), we double 
its capacity for doing work: 

Work = mx2hxg (1.1.2) 

Alternatively if we raise an object three times as heavy 
to a distance h above the ground (Figure 1.1. Id), the 
amount of work that this new object could perform would 
be three times that of the original object in Figure 1.1.1b: 

Work = 3mxhxg (1.1.3) 

The simple mechanical example in Figure 1.1.1 shows only 
one, simply understood way of doing work. Mechanical work 
can also be done by an object’s motion, as a demolition 
crew’s ‘wrecking ball’ illustrates. Electric current heating the 
element of an electric fire represents another form of work, 
as does an explosive charge used to blast a rock face in a 
mine. 



molecules in liquid water are more mobile than those in 
ice - that is, they have higher kinetic energy - the 
enthalpy of water (H wate ) is greater than that of an 
equivalent amount of ice (H. ce ) at the same temperature. 
The difference can be written: 

A H = H water -H lce (1.2) 

The A symbol (the Greek capital letter 'delta'), when 
written in front of H, signifies the difference in enthalpy 



between the initial (solid) and final (liquid) states of 
the compound the H,0. It represents the work (Box 1.1) 
that must be done in disrupting the chemical bonds 
that hold the crystal together. AH symbolizes the 
amount of heat that must be supplied from the sur- 
roundings for the crystal to melt completely; this is 
called the latent heat of fusion, or more correctly the 
enthalpy of fusion, a quantity that can be measured 
experimentally or looked up in tables. 
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(a) (b) (c) (d) 



Energy is simply the term that we use to describe a 
system’s capacity for doing work. Just as we recognize dif- 
ferent forms of work (mechanical, electrical, chemical ...), 
so energy exists in a number of alternative forms, as will 
be illustrated in the following pages. The energy stored in 
an electrical battery, for example, represents the amount 
of work that it can generate before becoming exhausted. 
A system’s capacity for doing work is necessarily expressed 
in the units of work (just as the capacity of a bucket is 
expressed as the number of litres of water it can contain), 
so it follows that energy is also expressed in 
joules = kg m 2 s’ 1 . When discussing large amounts of 
energy, we use larger units such as kilojoules (kJ = 10 3 J) 
or megajoules (MJ = 10 6 J). 



Figure 1.1.1 Work done in raising an object: (a) an object of 
mass m resting on the ground; (b) the same object elevated 
to height h; (c) the object elevated to height 2h; (d) another 
object of mass 3 m elevated to height h. Note-, elevation is 
measured between each object’s centre of gravity in its initial 
and final positions (note the centre of gravity of the larger 
weight is slightly higher than the smaller one). 



This simple example illustrates how one can go 
about documenting the energy changes that accom- 
pany geological reactions and processes, as a means of 
understanding why and when those reactions occur. 
This is the purpose of thermodynamics, a science that 
documents and explains quantitatively the energy 
changes in natural processes, just as economics anal- 
yses the exchange of money in international trade. 
Thermodynamics provides a fundamental theoretical 



framework for documenting and interpreting energy 
changes in processes of all kinds, not only in geology 
but in a host of other scientific disciplines ranging from 
chemical engineering to cosmology. 

Thermodynamics, because it deals with very abstract 
concepts, has acquired an aura of impenetrability in 
the eyes of many Earth science students, particularly 
those less at home in the realm of mathematics. With 
this in mind, one objective of these opening chapters 
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CHAPTER 1 



will be to show that thermodynamics, even at a simple 
and approachable level, can contribute a lot to our 
understanding of chemical reactions and equilibrium in 
the geological world. 

Energy changes in chemical systems are most easily 
introduced (as in Box 1.1) through analogy with 
mechanical forms of energy, which should be familiar 
from school physics. 

Energy in mechanical systems 



The energy of a body is defined as its capacity for doing 
work (Box 1.1). As we have discovered, 'work' can take 
various forms, but in simple mechanical systems it 
usually implies the movement of a body from one pos- 
ition to another against some form of physical resist- 
ance (friction, gravity, electrostatic forces, etc.). Then: 



Secondly, an object in a gravitational field possesses 
energy (i.e. can do work) by virtue of its position in 
that field, a property known as potential energy. The 
water held behind a hydroelectric dam has a high 
potential energy: under the influence of the Earth's 
gravitational field it would normally flow downhill 
until it reached sea-level, but the dam prevents this 
from happening. The fact that the controlled downward 
flow of this water under gravity can be made to drive 
turbines and generate electricity indicates that the 
water held behind the dam has the potential for doing 
work, and therefore possesses energy. 

The potential energy E p of an object of mass m at a 
height h above the ground is given by: 

E p = m x g x h (1 5) 

(J) ( k g)( m s~ 2 )( m ) 



work done = force required to move body 

J (joules) N ( nevvtons)=kg m s “ ^ 

x distance body is moved 

m (metres) 

So, for example, the work done in transporting a train- 
load of iron ore from A to B is the mechanical force 
required to keep the train rolling multiplied by the dis- 
tance by rail from A to B. The energy required to do this 
is generated by the combustion of fuel in the engine. 

One can distinguish two kinds of mechanical energy. 
Firstly, an object can do work by means of its motion. A 
simple example is the use of a hammer to drive a nail 
into wood. Work is involved because the wood resists 
the penetration of the nail. The energy is provided by 
the downward-moving hammer-head which, because 
of its motion, possesses kinetic energy (a name derived 
from the Greek kinetikos, meaning 'setting in motion'). 
The kinetic energy E k possessed by a body of mass m 
travelling with velocity v is given by: 



(J) ( k g)( m s _1 )‘ 

Thus the heavier the hammer (m) and/or the faster it 
moves (v), the more kinetic energy it possesses and the 
further it will drive in the nail. For similar reasons, a 
fast-moving stream can carry a greater load of sediment 
than a slow-moving one. 



where g is the acceleration due to gravity (9.81ms -2 ). 
Similar equations can be written representing the 
potential energies of bodies in other types of force 
field, such as those in electric and nuclear fields. 

An important aspect of potential energy is that its 
value as calculated from Equation 1.5 depends upon the 
baseline chosen for the measurement of height h. The 
potential energy calculated for an object in a second- 
floor laboratory, for example, will differ according to 
whether its height is measured from the laboratory floor, 
from the ground level outside, or from sea-level. The last 
of these alternatives seems at first sight to be the most 
widely applicable standard to adopt, but even that refer- 
ence point fails to provide a baseline that can be used for 
the measurement of height and potential energy down a 
deep mine (where both quantities might have negative 
values measured relative to sea-level). This ambiguity 
forces us to recognize that potential energy is not some- 
thing we can express on an absolute scale having a 
universal zero-point, as we do in the case of temperature 
or electric charge. The value depends upon the 'frame of 
reference' we choose to adopt. We shall discover that 
this characteristic applies to chemical energy as well. It 
seldom presents practical difficulties because in thermo- 
dynamics one is concerned with energy changes, from 
which the baseline-dependent element cancels out (pro- 
vided that the energy values used have been chosen to 
relate to the same frame of reference). 

In general, a body possesses kinetic and potential 
energy by virtue of its overall motion and position. 
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There is also an internal contribution to its total 
energy from the individual motions of its constituent 
atoms and molecules, which are continually vibrat- 
ing, rotating and - in liquids and gases - migrating 
about. This internal component, the aggregate of 
the kinetic energies of all the atoms and molecules 
present, is what we mean by the enthalpy of the body 
Enthalpy is closely related to the concept of heat (and 
was at one time referred to, rather misleadingly, as 
'heat content'). Heat is one of the mechanisms by 
which enthalpy can be transferred from one body to 
another. The effect of heating a body is simply to 
increase the kinetic energy of the constituent atoms 
and molecules, and therefore to increase the enthalpy 
of the body as a whole. 

Natural processes continually convert energy from 
one form into another. One of the fundamental axioms 
of thermodynamics, known as the First Law, is that 
energy can never be created, destroyed or 'lost' in such pro- 
cesses, but merely changes its form (Box 1.2). Thus the 
energy given out by a reaction is matched exactly by 
the amount of energy gained by the surroundings. 



Energy in chemical and mineral systems: 
free energy 



Experience tells us that mechanical systems in the 
everyday world tend to evolve in the direction that 
leads to a net reduction in total potential energy. Water 
runs downhill, electrons are drawn toward atomic 
nuclei, electric current flows from 'live' to 'neutral', 
and so on. The potential energy released by such 
changes reappears as other forms of energy or work: 
for example, the kinetic energy of running water, the 
light energy radiated by electronic transitions in atoms 
(Chapter 6), or the heat generated by an electric fire. 

Thermodynamics visualizes chemical processes in a 
similar way. Reactions in chemical or geological sys- 
tems arise from differences in what is called free energy, 
G, between products and reactants. The significance of 
free energy in chemical systems can be compared to 
that of potential energy in mechanical systems. A chem- 
ical reaction proceeds in the direction which leads to a 
net reduction in free energy, and the chemical energy so 
released reappears as energy in some other form - the 



' 

Box 1.2 The First Law of Thermodynamics 

The most fundamental principle of thermodynamics is that 
energy is never created, lost or destroyed. It can be trans- 
mitted from one body to another, or one place to another, 
and it can change its identity between a number of forms 
(as for example when the potential energy of a falling body 
is transformed into kinetic energy, or when a wind turbine 
converts the kinetic energy of moving air into electrical 
energy). But we never observe new energy being created 
from scratch, nor does it ever just disappear. Accurate 
energy bookkeeping will always show that in all known pro- 
cesses total energy Is always conserved. This cardinal prin- 
ciple is called the First Law of Thermodynamics. The energy 
given out by a reaction or process is matched exactly by 
the amount of energy gained by the surroundings. 

Implicit in the First Law is the recognition that work is 
equivalent to energy, and must be accounted for in energy 
calculations. When a compressed gas at room tempera- 
ture escapes from a cylinder, it undergoes a pronounced 
cooling, often to the extent that frost forms around the 
valve. (A smaller cooling effect occurs when you blow on 



your hand through pursed lips.) The cause of the cooling 
is that the gas has had to do work during escaping: it 
occupies more space outside the cylinder than when com- 
pressed inside it, and it must make room for itself by dis- 
placing the surrounding atmosphere. Displacing something 
against a resisting force (in this case atmospheric pres- 
sure) constitutes work, which the gas can only accomplish 
at the expense of its enthalpy. This is related directly to 
temperature, so that a drain on the gas’s internal energy 
reserves becomes apparent as a fall in temperature. 

A similar cooling effect may operate when certain gas- 
rich magmas reach high crustal levels or erupt at the sur- 
face. An example is kimberlite, a type of magma that 
commonly carries diamonds up from the mantle. Kimberlite 
penetrates to the surface from depths where the assoc- 
iated gases are greatly compressed, and the work that 
they do in expanding as the magma-gas system bursts 
through to the surface reduces its temperature; kimber- 
lites found in subvolcanic pipes (diatremes) appear to 
have been emplaced in a relatively cool state. 
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r 



Vapour 



v 







Water 



V J 

Figure 1.1 A simple model of chemical equilibrium 
between two coexisting phases, water and its vapour. 

The equilibrium can be symbolized by a simple equation. 

h 2 o ^ h 2 o 

liquid vapour 

At equilibrium, the migration of water molecules from the 
liquid to the vapour (evaporation, upward arrow) is 
balanced exactly by the condensation of molecules from 
vapour to liquid (downward arrow). 

electrical energy obtained from a battery, the light and 
heat emitted by burning wood, and so on. 

What form does free energy take? How can it be cal- 
culated and used? These questions are best tackled 
through a simple example. Imagine a sealed container 
partly filled with water (Figure 1.1). The space not 
filled by the liquid takes up water vapour (a gas) until 
a certain pressure of vapour is achieved, called the 
equilibrium vapour pressure of water, which is depend- 
ent only upon the temperature (assumed to be con- 
stant). H 2 0 is now present in two stable forms, each 
one having physical properties and structure distinct 
from the other: these two states of matter are called 
phases. From this moment on, unless circumstances 
change, the system will maintain a constant state, 
called equilibrium, in which the rate of evaporation 
from the liquid phase is matched exactly by the rate of 
condensation from the vapour phase: the relative 
volumes of the two phases will therefore remain 
constant. 

In this state of equilibrium, the free energies associ- 
ated with a given amount of water in each of these two 
phases must be equal. If that were not the case, a net 



flow of water molecules would begin from the high-G 
phase to the low-G phase, allowing the total free energy 
of the system to fall in keeping with the general ten- 
dency of chemical systems to minimize free energy. 
Any such flow, which would alter the relative propor- 
tions of the two phases, is inconsistent with the steady 
state observed. Clearly, at equilibrium, equivalent 
amounts of the two phases must have identical free 
energies: 



G 



vapour 



= G 



liquid 



( 1 . 6 ) 



This statement is in fact the thermodynamic definition 
of 'equilibrium' in such a system. 

But here we seem to have stumbled upon a paradox. 
Common sense tells us that to turn liquid water into 
vapour we have to supply energy, in the form of heat. 
The amount required is called the latent heat of evap- 
oration (more correctly, the enthalpy of evaporation). 
This indicates that the vapour has a greater enthalpy 
( H ) than an equivalent amount of the liquid 

v vapour' i i 

(^liquid) * 

H vapour > H l lq u ld O' 7 ) 



The difference reflects the fact that water molecules in 
the vapour state have (a) greater potential energy, hav- 
ing escaped from the intermolecular forces that hold 
liquid water together, and (b) greater kinetic energy 
(owing to the much greater mobility of molecules in 
the gaseous state). 

How can we reconcile Equations 1.6 and 1.7? Is it not 
common sense to expect the liquid state, in which the 
water molecules have much lower energies, to be 
intrinsically more stable than the vapour? What is it 
that sustains the vapour, in spite of its higher enthalpy, 
as a stable phase in equilibrium with the liquid? 

The answer lies in the highly disordered state char- 
acteristic of the vapour. Molecules in the gas phase fly 
around in random directions, occasionally colliding 
but, unlike molecules in a liquid, free to disperse 
throughout the available volume. The vapour is said to 
possess a high entropy (S). Entropy is a parameter that 
quantifies the degree of internal disorder of a substance 
(Box 1.3). Entropy has immense significance in thermo- 
dynamics, owing to Nature's adherence to the Second 
Law of Thermodynamics. This states that all spontaneous 
processes residt in an increase of entropy. The everyday 
consequences of the Second Law - so familiar that we 
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often take them for granted - are discussed further in 
Box 1.3. In the present context. Nature's preference for 
disordered, high-entropy states of matter is what 
makes it possible for vapour to coexist with liquid. In a 
sense, the higher entropy of the vapour 'stabilizes' it in 
relation to the liquid state, compensating for the higher 
enthalpy required to sustain it. 

Clearly, any analysis in energy terms, even of this 
simple example, will succeed only if the entropy differ- 
ence (AS) between liquid and vapour is taken into 
account. This is why the definition of the free energy 
(alternatively called the 'Gibbs energy') of each phase 
therefore incorporates an entropy term: 



^liquid 



— ^liquid f '^liquid 



( 1 . 8 ) 



C = Ft — T S 

'“vapour ^vapour “‘“vapour 



(1.9) 



H v and H are the enthalpies of the liquid and 

liquid vapour i i 

vapour respectively. S,. .. and S are the corres- 
ponding entropies. (Take care not to confuse the 
similar-sounding terms 'enthalpy' and 'entropy'.) The 
absolute temperature T (measured in kelvins) is 
assumed to be uniform in a system in equilibrium 
(Chapter 2), and therefore requires no subscript. 

The important feature of these equations is the nega- 
tive sign. It means that the vapour can have higher 
enthalpy ( H ) and higher entropy (S) than the liquid, 
and yet have the same free energy value (G), which 
must be true if the two phases are to be in equilibrium. 
Perhaps a more fundamental understanding of the 
minus sign can be gained by rearranging Equations 
1.8 and 1.9 into this form: 



Equations 1.8 and 1.9 express the fundamental contr- 
ibution that disorder makes to the energetics of chem- 
ical and geological reactions, a question we shall take 
up again in the following sections. 



Units 

Enthalpy, entropy and free energy, like mass and vol- 
ume, are classified as extensive properties. This means 
that their values depend on the amount of material 
present. On the other hand, temperature, density, vis- 
cosity, pressure and similar properties are said to be 
intensive properties, because their values are unrel- 
ated to the size of the system being considered. 

In published tables of enthalpy and entropy 
(Chapter 2), the values given are those for one mole - 
abbreviated in the SI system to 'mol' - of the substance 
concerned (18 g in the case of water). One therefore 
speaks of molar enthalpy and entropy, and of molar 
free energy and molar volume as well. The units of 
molar enthalpy and molar free energy are joules per 
mole (Jmol -1 ); those of molar entropy are joules per 
kelvin per mole (JK -1 mol -1 ). The most convenient 
units for expressing molar volume are 10 -6 m 3 mol -1 
(which are the same as cm 3 mol -1 , the units used in 
older literature). 

In thermodynamic equations like 1.8, temperature 
is always expressed in kelvins (K). One kelvin is equal 
in magnitude to one °C but the scale begins at the 
absolute zero of temperature (-273.15 °C), not at the 
freezing point of water (0 °C). Therefore: 



H = G + T.S (1.10) 

The enthalpy of a phase can thus be seen as consisting 
of two contributions: 



TlnK = Tln°C + 273.15 (1.11) 

The SI units for pressure are pascals (Pa; see 
Appendix A). 



G The part that potentially can be released by the 
operation of a chemical reaction, which is logic- 
ally called 'free' energy. This therefore provides 
a measure of the instability of a system (just as 
the potential energy of the water in a reservoir 
reflects its gravitational instability). 

T.S The part that is irretrievably bound up in the 
internal disorder of a phase at temperature T, 
and that is therefore not recoverable through 
chemical reactions. 



Free-energy changes 

For the reasons discussed above in relation to potential 
energy, the numerical values of G, . , and G have 

Oj" liquid vapour 

no absolute significance. In considering whether water 
will evaporate or vapour will condense in specific cir- 
cumstances, what concerns us is the change in free 
energy AG arising from the liquid-to-vapour 'reaction'. 
The first step in calculating free-energy changes is to 
write down the process concerned in the form of a 
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Box 1.3 Some properties of entropy 

The concept of disorder is of fundamental importance in thermodynamics, because it allows us to distinguish those pro- 
cesses and changes that occur naturally - ‘spontaneous’ processes - from those that do not. We are accustomed to 
seeing a cup shattering when it falls to the floor, but we never see the fragments reassemble themselves spontaneously 
to form a cup hanging on the dresser hook. Nor is it a natural experience for the air in a cold room to heat up a warm radia- 
tor. The direction of change that we accept as natural always leads to a more disordered state than we began with. 

To apply such reasoning to the direction of chemical change, we need a variable that quantifies the degree of 
disorder in a chemical system. In thermodynamics this is defined by the entropy of the system. To define entropy 
rigorously lies beyond the scope of this book, but it is worth identifying the processes that lead to an increase of 
entropy. The entropy of a system depends upon: 

(i) the distribution of matter or of individual chemical species in the system; and 

(ii) the distribution of energy. 



Entropy and the distribution of matter 

(a) Entropy increases as a substance passes from the solid state to the liquid state to the gaseous state.* 

(a) 



Crystal Liquid/glass* Gas 

(b) Entropy increases when a gas expands. 



(b) 



O o ° 0 
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Ordered state Disordered state 



(c) Entropy increases when pure substance are mixed together. 



(c) 




Ordered state 



Disordered state 



* A glass is a solid having the disordered structure of a liquid, but deprived of atomic mobility (no flow). Its entropy is intermediate 
between liquid and crystalline solid. 

V 
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Entropy and the distribution of energy 

The entropy of a system increases 

(d) 



when a substance is heated: 



because raising the temperature makes the random thermal motions of atoms 
and molecules more vigorous. 



(e) 

when heat flows from a hot body (e.g. a radiator) to a cold body (the surrounding air): 




Heat concentrated in one place 

Ordered state 



Heat uniformly distributed 

Disordered state 




(f) 

when chemical energy (fuel + oxidant) is transformed into heat: 





Chemical energy of Hot combustion products 

fuel/oxidant mixture Disordered state 

Ordered state 

(g) 

when mechanical enerav is transformed into heat (e.o. frictionl. 



The Second Law of Thermodynamics 

The Second Law states that the operation of any spontaneous process leads to an (overall) increase of entropy. Our experi- 
ence of this law is so intricately woven into the fabric of everyday life that we are scarcely aware of its existence, but its 
impact on science is nonetheless profound. The expansion of a gas is a spontaneous process involving an increase of 
entropy: a gas never spontaneously contracts into a smaller volume. Water never runs uphill. Applying heat to an electric 
fire will never generate electricity. All of these impossible events, were they to occur, would bring about a reduction of 
entropy and therefore violate the Second Law. 

Entropy is lowest when energy is concentrated in one part of a system. This is a characteristic of all of the energy resources 
that we exploit: water retained behind a hydroelectric dam, chemical energy stored in a tank of gasoline or in a charged battery, 
nuclear energy in a uranium fuel rod, etc. Entropy is highest when energy is evenly distributed throughout the system being 
considered, and in such circumstances it cannot be put to good use. Spontaneous (entropy-increasing) changes are always 
accompanied by a degradation in the ‘quality’ of energy, in the sense that it becomes dispersed more widely and uniformly. 









10 



CHAPTER 1 



chemical reaction, as in Equation 1.1. For the water/ 
vapour equilibrium: 

H 2 0 ^ H z O (1.12) 

liquid vapour 

The equilibrium symbol (^) represents a balance 
between two competing opposed 'reactions' taking 
place at the same time: 



'Forward' reaction: liquid vapour (evaporation) 
reactant product 

'Reverse' reaction: liquid vapour (condensation) 

product reactant 

By convention, the tree-energy change for the forward 
reaction (AG) is written: 



AG G products G reacdln j s 
= G — G 

'""vapour '"'liquid 



(1.13) 



Each G can be expressed in terms of molar enthalpy 
and entropy values obtained from published tables 
(Equations 1.8 and 1.9). Thus 



_T C 

vapour ‘‘"vapour 



AG = (H 

-(»■ 

= AH-T.AS 



) (^liquid "GSliquidj 



— H 

vapour ' liquid 



)-T(S 



- s 

vapour ‘"'liquid J 



(1.14) 



In this equation AH is the heat input per mole required 
to generate vapour from liquid (the latent heat of 
evaporation). In the context of a true chemical reac- 
tion, it would represent the heat of reaction (strictly the 
enthalpy of reaction). If AH for the forward reactions 
is negative, heat must be given out by the reaction, 
which is then said to be exothermic ('giving out 
heat'). A positive value implies that the reaction will 
proceed only if heat is drawn in from the surround- 
ings. Reactions that absorb heat in this way are said 
to be endothermic ('taking in heat'). AS represents 
the corresponding entropy change between liquid 
and vapour states. 

The values of H . H, ... S and S, . , can be 
looked up as molar quantities for the temperature of 
interest (e.g. room temperature ~ 298 K) in published 
tables. In this case, AH and AS can be calculated by 
simple difference, leading to a value for AG (taking 
care to enter the value of T in kelvins, not °C). From 
the sign obtained for AG, it is possible to predict in 



which direction the reaction will proceed under the 
conditions being considered. A negative value of AG 
indicates that the products are more stable - have a 
lower free energy - than the reactants, so that the 
reaction can be expected to proceed in the forward 
direction. If AG is positive, on the other hand, the 
'reactants' will be more stable than the 'products', 
and the reverse reaction will predominate. In either 
case, reaction will lead eventually to a condition 
where AG = 0, signifying that equilibrium has been 
reached. 

Now let us see how these principles apply to miner- 
als and rocks. 



Stable, unstable and metastable minerals 



The terms 'stable' and 'unstable' have a more precise 
connotation in thermodynamics than in everyday 
usage. In order to grasp their meaning in the context of 
minerals and rocks, it will be helpful to begin by con- 
sidering a simple physical analogue. Figure 1.2a shows 
a rectangular block of wood in a series of different 
positions relative to some reference surface, such as a 
table top upon which the block stands. These config- 
urations differ in their potential energy, represented by 
the vertical height of the centre of gravity of the block - 
shown as a dot - above the table top. Several general 
principles can be drawn from this physical system 
which will later help to illuminate some essentials of 
mineral equilibrium: 

(a) Within this frame of reference, configuration D has 
the lowest potential energy possible, and one calls 
this the stable position. At the other extreme, con- 
figurations A and C are evidently unstable, because 
in these positions the block will immediately fall 
over, ending up in a position like D. Both clearly 
have higher potential energy than D. 

(b) In discussing stable and unstable configurations, 
one need not consider all forms of energy pos- 
sessed by the wooden block, some of which 
(for example, the total electronic energy) would 
be difficult to quantify. Mechanical stability 
depends solely upon the relative potential energies 
of - or energy differences between - the several 
configurations, and not on their absolute energy 
values. 
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Figure 1.2 Potential energy of a rectangular wooden block in various positions on a planar surface, (a) Four positions of the 
block, showing the height of its centre of gravity (dot) in each case, (b) The pattern of potential-energy change as the block 
topples, for the unstable (A, C) and metastable (B) configurations. 



(c) Configuration B presents something of a paradox. 
It has a potential energy greater than the unstable 
state C, yet, if left undisturbed, it will persist indef- 
initely, maintaining the appearance of being stable. 
The introduction of a small amount of energy, such 
as a person bumping into the table, may however 
be sufficient to knock it over. The character of con- 
figuration B can be clarified by sketching a graph 
of potential energy against time as the block top- 
ples over (Figure 1.2b). For both unstable positions 
A and C, the potential energy falls continuously to 
the value of position D; but in the case of position 
B the potential energy must first rise slightly, before 
falling to the minimum value. The reason is that 
the block has to be raised on to its corner (position 
A) before it can fall over, and the work involved in 
so raising its centre of gravity constitutes a poten- 
tial energy 'hurdle' which has to be surmounted 
before the block can topple. By inhibiting the spon- 
taneous toppling of the block, this hurdle stabilizes 
configuration B. One uses the term metastable to 
describe any high-potential-energy state that is 
stabilized by such an energy hurdle. 



The application of this reasoning to mineral stability 
can be illustrated by the minerals calcite and aragonite, 
whose ranges of stability in pressure-temperature 
space are shown in the form of a phase diagram in 
Figure 1.3a. These minerals are alternative crystallo- 
graphic forms of calcium carbonate (CaC0 3 ), stable 
under different physical conditions. The phase diagram 
shown in Figure 1.3a is divided into two areas called 
stability fields, one representing the range of applied 
pressure and temperature under which calcite is the 
stable mineral; the other - at higher pressures - indi- 
cating the range of conditions favouring aragonite. 
The stability fields are separated by a line, called a 
phase boundary, which defines the restricted set of 
circumstances under which calcite and aragonite can 
coexist together in equilibrium with each other. 

The energetics of the calcite-aragonite system are 
illustrated in Figure 1.3b, which shows how the molar 
free energies of the two minerals vary along the line 
X-Y in Figure 1.3a. At high pressure (Y), deep within 
the crust, the molar free energy of aragonite is less than 
that of calcite, and thus aragonite is the stable mineral 
under these conditions, analogous to configuration D 
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(a) 




(b) 




Figure 1.3 Stability of CaC0 3 polymorphs, (a) Pressure-temperature phase diagram showing the stability fields of calcite 
and aragonite. The line A-B is the phase boundary, indicating the P-T conditions under which calcite and aragonite are able 
to coexist in stable equilibrium, (b) Variation of free energy (G) of calcite and aragonite over a range of pressures along the 
isothermal line X-Y in (a). 



of the wooden block in Figure 1.2a. At a lower pressure 
(X) nearer to the surface, however, the position is 
reversed: calcite has the lower free energy and is there- 
fore the stable mineral. The lines representing the free 
energy of calcite and aragonite as a function of pres- 
sure cross over in Figure 1.3b at a point marked p. Here 
the two minerals have equal molar free energies, and 
are therefore in chemical equilibrium with each other. 
Point p therefore marks the position in Figure 1.3b of 
the phase boundary appearing in Figure 1.3a. 

Imagine transporting a sample of aragonite from the 
conditions represented by point Y to a new location (at a 
shallower depth in the Earth's crust) having the pressure 
and temperature coordinates of point X. Under the new 
conditions aragonite will no longer be the stable mineral 
and it will tend to achieve a state of lower free energy by 
recrystallizing to calcite. This transformation may not 
occur immediately, however, because the status of arag- 
onite is similar to the wooden block in position B. The 
same three points made above in relation to Figure 1.2 
can be reiterated for the calcite-aragonite system: 

(a) Calcite, having the lower free energy, will be the 
stable form of calcium carbonate under the lower 
pressure conditions defined by X. 



(b) Many other forms of energy are associated with 
calcite and aragonite under such conditions, but in 
discussing thermodynamic stability we are con- 
cerned only with free-energy differences between 
alternative states. This has the important conse- 
quence that free energy needs only to be expressed 
in relative terms, referred to a convenient but 
arbitrary common point, a sort of thermodynamic 
'sea-level'. All important applications of thermo- 
dynamics involve the calculation of free-energy 
differences between the various states of the sys- 
tem being considered, and the notion of an absol- 
ute scale of free-energy values, analogous to the 
absolute temperature scale, is unnecessary and 
inappropriate. 

(c) In spite of not being stable under near-surface con- 
ditions (Figure 1.3), aragonite is quite a common 
mineral, and may survive for long periods of geol- 
ogical time on the surface of the Earth. Like con- 
figuration B in Figure 1.2a, aragonite may give the 
appearance of being in a stable state in such cir- 
cumstances, even though its free energy clearly 
exceeds that of calcite. (Under certain circum- 
stances, aragonite may actually crystallize under 
near-surface conditions: for example, the shells of 
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planktonic gastropods (pteropods) are constructed 
of aragonite precipitated directly from seawater.) 
The explanation for this appearance of stability is 
that aragonite is stabilized by a tree-energy 'hur- 
dle', just like its mechanical analogue. 

The tree-energy path followed by aragonite (under 
the conditions shown by X in Figure 1.3a), as it under- 
goes transformation into calcite, is shown Figure 1.4. 
The energy hurdle exists because rearranging the 
crystal structure of aragonite into that of calcite 
involves some work being done in breaking bonds 
and moving atoms about. Although this energy invest- 
ment is recovered several times over in the net release 
of free energy as the reaction proceeds, its importance 
in determining whether the reaction can begin is con- 
siderable. In recognition of this influence, the height 
of the energy hurdle is called the activation energy of 
the reaction (symbol E a ). 

Figure 1.4 illustrates an important distinction 
between two major domains of geochemical science. 
Thermodynamics is concerned with the tree-energy 
changes associated with chemical equilibrium between 
phases, and provides the tools for working out which 
mineral assemblages will be stable under which condi- 
tions. Only the initial and final states are of interest in 
thermodynamics, and attention is confined to net 
energy differences between reactants and products 
(AG, AH, AS), as Chapter 2 will show. The science of 
chemical kinetics deals with the mechanics of the reac- 
tions that lead to equilibrium and the rates at which 
they occur. In this area, as we shall see in Chapter 3, the 




Progress of reaction 



Figure 1.4 The tree-energy path during the recrystallization 
of aragonite into calcite at the P-T conditions shown by X in 
Figure 1.3. 

activation energy asserts a dominant role, accounting 
for the strong influence of temperature on many 
geological processes. 
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2 EQUILIBRIUM IN GEOLOGICAL 
SYSTEMS 



The significance of mineral stability 



Igneous and metamorphic rocks form in places that are 
not, on the whole, directly accessible to the investigat- 
ing geologist. To discover how such rocks are produced 
within the Earth, one must resort to indirect lines of 
enquiry. The most important clues consist of the min- 
erals that the rocks themselves contain. A given min- 
eral crystallizes as a stable phase only within restricted 
ranges of pressure and temperature, as we saw in the 
case of aragonite which occurs stably only at high pres- 
sures (Figure 1.3a). Subjecting the mineral to condi- 
tions that fall outside its stability range will cause 
another mineral, which is stable under those condi- 
tions (like calcite at low pressure), to begin crystalliz- 
ing in its place. The stability of other types of minerals 
may depend in a similar way on the pressure of water 
vapour or some other gaseous component present 
during crystallization, and such a mineral will occur in 



a rock only if the vapour pressure present during its 
formation falls within the appropriate range. 

The sensitivity of such minerals to the physical cir- 
cumstances of their formation offers the petrologist 
tremendous opportunities because, when found in an 
igneous or metamorphic rock now exposed at the sur- 
face, they provide a means of establishing quantita- 
tively the characteristics of the physical environment 
in which that rock originally crystallized. The study of 
mineral stability therefore offers the key to a veritable 
library of information, sitting in the rocks waiting to be 
utilized, about conditions and processes deep within 
the Earth's crust and upper mantle. 

The usual way to establish the physical limits within 
which a mineral is stable - and beyond which it is 
unstable - is to cook it up in a laboratory experiment 
(Box 2.1). Technology today is capable of reproducing in 
the laboratory the physical conditions (temperature 
T, pressure P, water vapour pressure P H Q , and so on) 



Chemical Fundamentals of Geology and Environmental Geoscience, Third Edition. Robin Gill. 
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Box 2.1 Phase equilibrium experiments with minerals 



The process of mapping out a phase diagram like that in 
Figure 2.1 is illustrated in Figure 2.1.1a. Each of the filled 
and open circles represents an individual experiment in 
which a sample of the relevant composition is heated in a 
pressure-vessel at the pressure and temperature indicated 
by the coordinates, for a sufficient time for the phases to 
react and equilibrate with each other. At the end of each 
experiment - which may last hours, days or even months, 
depending upon the time needed to reach equilibrium - the 
sample is ‘quenched’, meaning that it is cooled as quickly 
as possible to room temperature in order to preserve the 
phase assemblage formed under the conditions of the 
experiment (which on slower cooling might recrystallize to 
other phases - see Chapter 3). The sample is removed 
from its capsule, and the phase assemblage is identified 
under the microscope or by other methods. The symbol for 
each experiment is ornamented on the diagram in such a 
way that it indicates the nature of the phase assemblage 
observed, so that the results of a series of experiments 
allow the position of the phase boundary to be determined. 
Conditions can be chosen for later experiments which 
allow accurate bracketing of its position in P-T space. 



Experiments in the laboratory must necessarily be 
concluded in much shorter times than nature can take to 
do the same job. Even at high temperatures, silicate 
reactions are notoriously sluggish, and the assemblage 
observed at the end of an experimental run might reflect 
an incomplete reaction or a metastable intermediate 
state rather than a true equilibrium assemblage. The pro- 
portions that this problem can sometimes assume are 
illustrated by the disagreement among the published 
determinations of the kyanite-sillimanite-andalusite triple 
point shown in Figure 2.1.1b. The present consensus 
places the triple point at about 4 x 10 s Pa and 500 °C 
(Figure 2.1). 

One precaution that the experimenter can take is to 
ensure that the position of every phase boundary is estab- 
lished by approaching it from both sides, a procedure 
known as ‘reversing the reaction’. In locating the kyanite- 
sillimanite phase boundary, for example, it is insufficient 
just to measure the temperature at which kyanite changes 
into sillimanite; the careful experimenter will also meas- 
ure the temperature at which sillimanite, on cooling, 
inverts to kyanite. 
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Figure 2.1.1 (a) A P-T phase diagram showing the reaction curve for reaction 2.5. Each dot represents the P-T values 
for an individual experiment. All experiments were conducted in a C0 2 atmosphere; the pressure of C0 2 gas present 
(symbolized as P C02 ) is equal to the applied pressure. Solid dots identify runs that produced calcite and quartz; 
open dots represent runs in which wollastonite was formed. The curve is drawn to run between filled and open dots. 

(b) Differences between published P-T values of the kyanite-sillimanite-andalusite triple point shown in Figure 2.1. 
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encountered anywhere in the crust and in most of the 
upper mantle. Experiments can be conducted in which 
minerals are synthesized at a series of accurately known 
temperatures and pressures, from starting materials of 
known composition. In each case, the phases (minerals) 
formed are considered to crystallize in chemical equilib- 
rium with each other, and/or with a molten silicate 
liquid, and for this reason the experiments are called 
phase-equilibrium experiments. Experimental information 
like this is vital to the petrological interpretation of 
naturally occurring rocks. The results are used to map 
out, in phase diagrams similar to Figure 1.3a, the areas 
within which particular mineral assemblages are stable. 
In crossing a boundary from one stability field into a 
neighbouring one as physical conditions change, one 
mineral assemblage recrystallizes into another, a chem- 
ical reaction between coexisting minerals which trans- 
forms an unstable - or metastable - assemblage into a 
new, stable one. These boundaries, like the diagonal 
line in Figure 1.3a, are called phase boundaries or, more 
generally, reaction boundaries. 

Phase diagrams form a key part of the literature of 
petrology, a powerful means of portraying and inter- 
preting the phase-equilibrium data relevant to igneous 
and metamorphic petrogenesis. Reading and inter- 
preting such diagrams, in the light of underlying 
thermodynamic principles, is one of the basic skills 
essential to every geologist. 

Just as atmospheric pressure represents the weight 
of the column of air above us under the influence the 
Earth's gravitational field, so the pressure experienced 
by a rock buried at some depth within the Earth reflects 
the total weight of the column of rock + ocean + atmos- 
phere resting upon it. This lithostatic pressure P there- 
fore increases with depth d in a simple manner 
conveniently approximated as: 

P ~3xd (2.1) 

where P is expressed in units of 10 8 Pa and d is in km. 
The high pressure applied in some phase-equilibrium 
experiments is simply the means by which we simu- 
late in the laboratory the effect of depth in the Earth. 

Systems, phases and components 



To avoid confusion, one must be clear about the mean- 
ing of several terms which are used in a specific sense 
in the context of phase equilibrium. 



System 

This is a handy word for describing any part of the 
world to which we wish to confine attention. 
Depending on the context, system could mean the 
whole of the Earth's crust, or the oceans, or a cooling 
magma chamber, or an individual rock, or a sample 
undergoing a phase-equilibrium experiment. In most 
cases the term will refer to a collection of geological 
phases (see below) interacting with each other. 

An open system is one that is free to exchange both 
material and energy with its surroundings. The sea is 
an open system; it may be useful to consider it as an 
isolated entity for the purposes of discussion, but one 
must recognize that it receives both sunlight and river 
water from outside, and loses heat, water vapour and 
sediment to the atmosphere and crust. 

A closed system is one that is sealed with respect to 
the transfer of matter, but that can still exchange energy 
with the surroundings. A sealed magma chamber 
would be a good example, its only interaction with its 
environment (ideally) being the gradual loss of heat. 
An isolated system is one that exchanges neither mass 
nor energy with its surroundings, a notion of little 
relevance to the real world. 

'System' may alternatively be used to denote a domain 
of chemical (rather than physical) space. Petrologists use 
the word to distinguish a particular region of 'composi- 
tional space' to which attention is to be confined. Thus 
one speaks of 'the system MgO-SiO z ', referring to the 
series of compositions that can be generated by mixing 
these two chemical components in all possible pro- 
portions. The system so designated includes various 
minerals whose compositions lie within this range (the 
silica minerals, and olivine and pyroxene). 

Phase 

The meaning of phase in physical chemistry and petrol- 
ogy is easy to grasp although cumbersome to put into 
words. In formal terms, a phase can be defined as 'a 
part or parts of a system occupying a specific volume 
and having uniform physical and chemical characteristics 
which distinguish it from all other parts of the system'. 

Each individual mineral in a rock thus constitutes a 
separate phase, but that is not the end of the story. A 
lump of basalt collected from a solidified lava flow 
might be found to consist of four minerals, say plagio- 
clase, augite, olivine and magnetite. But the igneous 
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texture tells us that these minerals crystallized from a 
melt, which itself had 'uniform physical and chemical 
properties' distinct from each of the crystalline miner- 
als present. So in considering the phase relations that 
dictated the present character of the rock we must 
count the melt as a phase too, although it is now no 
longer present as a constituent. If the basalt is vesicular 
(that is, it contains empty bubbles or 'vesicles'), we 
have evidence that a sixth phase, water vapour, was 
also present as the rock crystallized. 

(Note that water can be present as a dissolved spe- 
cies within a melt or a hydrous crystalline mineral. 
Being accommodated within the volume occupied by 
that phase, it does not then count as a phase in its own 
right. Only when separate bubbles of water vapour 
appear can one accord water the status of a separate 
phase. When this happens, as with the vesiculating 
lava, it is a sign that all other phases present contain as 
much water as they can accommodate, and the system 
is saturated with water.) 

One must therefore be careful not to overlook addit- 
ional phases which - though no longer present in a 
particular rock - may have influenced its formation or 
development. Here are some examples where this 
could happen: 

(a) Some igneous rocks show textural evidence for the 
existence of two distinct (presumably immiscible) 
silicate liquids that once existed together in mutual 
equilibrium. 

(b) A metamorphic rock may have developed in the 
presence of a vapour phase, permeating the grain 
boundaries between crystals, of which no visible 
trace now remains. 

(c) Mineral veins are deposited from a fluid phase 
preserved only in occasional microscopic inclu- 
sions within crystals (Box 4.6). 

(d) When a magma is produced by melting deep 
inside the Earth, it may be in pressure-dependent 
equilibrium with minerals quite different from 
those that crystallize from the same melt at the 
surface. 

Thus the minerals we see in a rock today may repre- 
sent only a part (or even no part at all) of the original 
phase equilibria to which the rock owes its present 
constitution. 

It is usual to refer to solid phases by the appropriate 
mineral name - quartz, kyanite, olivine, and so on. The 



essential distinction between them is crystallographic 
structure, not their chemical composition (which sev- 
eral minerals may have in common). Any molten phase 
present - regardless of composition - is referred to as 
'melt'. By convention one refers to any gaseous phase 
involved as 'vapour'. 

Component 

The basic chemical constituents of a system, of which 
the various phases are composed, are called its compo- 
nents. The concept of a component is defined in a pre- 
cise but rather roundabout way: 'the components of a 
system comprise the minimum number of chemical 
(atomic and molecular) species required to specify 
completely the compositions of all the phases 
present'. 

Consider a crystal of olivine, which at its simplest 
consists of the elements magnesium (chemical symbol 
Mg), iron (Fe), silicon (Si) and oxygen (O). One way to 
define the components of the olivine would be to 
regard each chemical element as a separate component 
because the composition of any olivine can be stated in 
terms of the concentrations 1 of four elements: 

Mg Fe Si O 

However, defining the components in this way 
fails to recognize an important property of all sili- 
cate minerals, including olivine: that the oxygen 
content is not an independent quantity, but is tied by 
valency (Chapter 6) to the amounts of Mg, Fe and Si 
present, being just sufficient to generate the oxides 
of each of these elements (this is explained in 
Box 8.4). So, in describing the composition of an 
olivine the same information can be conveyed more 
economically in terms of the concentrations of only 
three components: 

MgO FeO SiO z 

By using a property specific to olivines, however, a 
still more economical statement of olivine composition 
can be devised. The crystal chemistry of olivine (see 
Chapter 8) requires an olivine composition to conform 
to a general formula that we can represent by X,SiO ( . X 
represents a type of atomic site in the olivine crystal 

1 The ways in which concentration can be expressed are 

summarized at the beginning of Chapter 4. 
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structure that accommodates either Mg or Fe, but not 
Si. For every Si atom in the olivine structure, there have 
to be two divalent atoms present, each of which can be 
either Mg or Fe. Another way to symbolize this con- 
straint is to write the formula as (Mg,Fe) 2 Si0 4 , in which 
'(Mg,Fe)' represents an atom of either Mg or Fe. One 
can now express the composition of an olivine as a 
combination of just two components: 

Mg 2 Si0 4 Fe 2 Si0 4 

Mineralogists call these components the 'end-mem- 
bers' of the olivine series, and give them the names 
forsterite (Fo) and fayalite (Fa) respectively. 

In analysing the arithmetic of chemical equilib- 
rium between minerals, it is important to formulate 
the components of a system in such a way as to mini- 
mize their number, as the definition implies. What 
constitutes the minimum number depends upon the 
nature of the system. In an experiment involving 
the melting of an olivine crystal on its own, the comp- 
osition of the melt, though different from the solid, still 
conforms to the olivine formula X,SiC) 4 . The comp- 
ositions of both phases present, olivine and melt, can 
therefore be expressed as proportions of only two 
components, Mg 2 Si0 4 and Fe,SiC) 4 (Box 2.4). Systems 
consisting of only two components are called binary 
systems. 

If, however, olivine coexists with, let us say, orthopyr- 
oxene, the formulation of components becomes less 
straightforward. Orthopyroxenes are composed of the 
same four elements as olivine, but they combine in 
different proportions. The general formula of ortho- 
pyroxene, X 2 Si 2 0 6 , reveals an X:Si ratio (1:1) lower than 
for olivine (2:1). The composition of a pyroxene cannot 
therefore be expressed in terms of just the two olivine 
end-members. To represent the separate compositions 
of olivine and orthopyroxene in this system, three 
components will be needed: 

Either : Mg 2 Si0 4 Fe 2 Si0 4 SiO z 
Or : MgO FeO SiO z 

The identity of the components is less important here 
than their number. A system like this requiring three 
components to express all possible compositions is 
said to be ternary. 

There are circumstances in which four components 
would be necessary, such as when olivine coexists with 



metallic iron (for example in certain meteorites). The 
amount of oxygen present is no longer determined 
solely by the metals present, as it would be in a system 
consisting entirely of silicates. One cannot express the 
composition of metallic iron as a mixture of oxides, so 
one must resort to using four components. Mg, Fe, Si 
and O, in order to describe all possible compositions in 
this quaternary system. 

In this book, the general practice will be to refer to 
components by means of their chemical formulae. This 
avoids confusion between phases and components, 
which can arise when a phase (for example the mineral 
quartz) happens to have the same chemical compos- 
ition as one of the components (SiO,) in the same sys- 
tem. However, in other books it is quite common for 
end-member names to be used in this way as well (for 
example, 'forsterite' for Mg,Si0 4 ). 

Equilibrium 



It is useful to distinguish between two aspects 
of equilibrium: thermal equilibrium and chemical 
equilibrium. 

Thermal equilibrium 

All parts of a system in thermal equilibium have the 
same temperature: in these circumstances heat flowing 
from one part of the system. A, to another part, B, is 
exactly balanced by the heat passing from part B to 
part A, so there is no net transfer of heat. Net heat 
transfer only occurs when there is a difference in tem- 
perature between different parts of the system. 

Chemical equilibrium 

This describes a system in which the distribution of 
chemical components among the phases of a system 
has become constant, showing no net change with 
time. This steady state does not mean that the flow of 
components from one phase to another has ceased: 
equilibrium is a dynamic process. An olivine sus- 
pended in a magma is constantly exchanging comp- 
onents with the melt. At melt temperatures, atoms will 
diffuse across the crystal boundary, both into the crys- 
tal and out of it into the liquid. If the diffusion rates of 
element X in and out of the crystal are unequal, there 
will be a net change of the composition of each phase 
with time, a condition known as disequilibrium. Such 
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changes usually lead eventually to a situation where, 
for every element present, the flux of atoms across the 
crystal boundary is the same in both directions, result- 
ing in zero net flow, and no change of composition 
with time. This is what we mean by equilibrium. 

The rate at which equilibrium is achieved varies 
widely and, as Chapter 3 will show, disequilibrium is 
found to be a common condition in geological systems, 
particularly at low temperatures. 

The Gibbs Phase Rule 

A natural question to ask is: how many phases can be 
in equilibrium with each other at any one time? In 
Figure 1.3a we looked at a simple system in which only 
two phases occurred. Most actual rocks, however, are 
not so simple. What factors determine the mineralog- 
ical complexity of a natural rock? Which aspect of a 
chemical equilibrium controls the number of phases 
that participate in it? 

This question was addressed in the 1870s by the 
American engineer J. Willard Gibbs, the pioneer of 
modem thermodynamics. The outcome of his work 
was a simple but profoundly important formula called 
the Phase Rule, which expresses the number of phases 
that can coexist in mutual equilibrium (</>) in terms of 
the number of components (C) in the system and 
another property of the equilibrium called the variance 
(F). The Phase Rule can be stated symbolically as: 

f + F = C + 2 (2-2) 

The variance is alternatively known as the number of 
degrees of freedom (hence the symbol F used to represent 
it). The concept is most easily introduced through an 
example. Figure 2.1 illustrates the equilibrium phase 
relations between the minerals kyanite, sillimanite and 
andalusite. These minerals are all aluminium silicate 2 
polymorphs of identical composition. A single compo- 
nent (Al 2 SiO s ) is therefore sufficient to cover the com- 
positional 'range' of the entire system. 

Points A, B and C are three different points in the 
'P-T space' covered by the diagram; they represent 
three classes of equilibrium that can develop in the sys- 
tem. The obvious difference between them is the nature 



2 Not to be confused with aluminosilicate minerals discussed 
in Chapter 9. 




Figure 2.1 A P-T diagram showing phase relations between 
the aluminium silicate minerals (composition Al 2 SiO s ). The 
pressure axis is graduated in units of 10 s pascals, equal in 
magnitude to the traditional pressure units, kilobars 
(lkbar = 10 3 bars = 10 8 Pa). Kyanite is a triclinic mineral that 
usually occurs as pale blue blades in hand specimen. 
Sillimanite is orthorhombic and is commonly fibrous or 
prismatic in habit. Andalusite is also orthorhombic and is 
characteristically pink in hand specimen. 

of the equilibrium assemblage. Point A lies within a 
field where only one phase, kyanite, is stable. Point B 
lies on the phase boundary between two stability fields, 
where two minerals, kyanite and sillimanite, are stable 
together. Point C, at the triple point where the three 
stability fields (and the three phase boundaries) meet, 
represents the only combination of pressure and tem- 
perature in this system at which all three phases can 
exist stably together. 

It is clear that the three-phase assemblage 
(kyanite + sillimanite + andalusite), when it occurs, 
indicates very precisely the state of the system (that is, 
the values of P and T) in which it is produced, because 
there is only one set of conditions under which this 
assemblage will crystallize in equilibrium. Using the 
Phase Rule (Equation 2.2), one can calculate that the 
variance F at point C is zero: 



Point C cj)=3 


(3 phases, ky+sili + andal) 


C= 1 

3+F=l+2 


(1 component, AI 2 Si0 5 ) 


Therefore F= 0 


an invariant equilibrium 
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A variance of zero means that the three-phase equilib- 
rium assemblage completely constrains the state of the 
system to a particular combination of P and T. Such a 
situation is called invariant. There is no latitude (no 
degree of freedom) for P or T to change at all, if the 
assemblage is to remain stable. Any such variation 
would lead to the disappearance of one or more of the 
three phases, thus altering the character of the equilib- 
rium. So the discovery of this mineral assemblage in a 
metamorphic rock (a rare event, since only about a 
dozen natural occurrences are known) ties down very 
precisely the conditions under which the rock must 
have been formed, assuming that: 

(a) the kyanite-andalusite-sillimanite assemblage rep- 
resents an equilibrium state obtained as the rock 
formed, and not simply an uncompleted reaction 
from one assemblage to another; and 

(b) the P-T coordinates of the invariant point in 
Figure 2.1 are accurately known from experimental 
studies. (Whether this requirement is satisfied in 
the case of the Al 2 SiO s polymorphs is arguable - 
see Box 2.1 - but in the present discussion we shall 
ignore this difficulty.) 

The two-phase equilibrium between, say, kyanite 
and sillimanite is less informative. The coexistence of 
these two minerals indicates that the state of the system 
in which they crystallized must lie somewhere on the 
kyanite-sillimanite phase boundary, but exactly where 
along this line remains uncertain unless we can specify 
one of the coordinates of point B. We only need to spec- 
ify one coordinate (for example, temperature) because 
the other is then fixed by intersection of the specified 
coordinate with the phase boundary. According to the 
Phase Rule, the variance at B is equal to 1: 



Point B (f> = 2 


(2 phases, ky+sill) 


C=1 


(1 component, Al 2 Si0 5 ) 


2+F=l+2 




Therefore F= 1 


a univariant equilibrium 



The one degree of freedom indicates that the state of the 
system is only unconstrained in one direction in 
Figure 2.1, along the phase boundary. One additional 
piece of information is required (either T or P) to tie 
down the state of the system completely. The coexistence 
of kyanite and sillimanite in a rock will pinpoint the 



exact conditions of origin only in conjunction with other 
information about P or T. 

At point A, where kyanite occurs alone, the variance 
is equal to 2: 



Point A <j> = l 


(1 phase, ky) 


C= 1 

l+F=l+2 


(1 component, AI 2 Si0 5 ) 


Therefore F= 2 


a divariant condition 



Within the bounds of the divariant kyanite field, there- 
fore, P and T can vary independently (two degrees of 
freedom) without upsetting the equilibrium phase 
assemblage (just kyanite). The one-phase assemblage 
is therefore little help in establishing the precise state 
of the system, because it leaves two variables (P and T) 
still to be specified. 

The variance cannot be greater than 2 in a one-com- 
ponent system like Figure 2.1. In the more complex 
systems we shall meet in the following section, the 
phases present may consist of different proportions of 
several components. A complete definition of the state 
of such a system must then include the compositions of 
one or more phases, in addition to values of P and T. 
Such compositional terms (X a , X h , etc., representing the 
mole fractions of a, b, etc. in a phase) contribute to the 
total variance, which can therefore, in multicomponent 
systems, adopt values greater than 2. 

Variance can be summarized in the following way. 
The 'state' of a system - whether we consider a simple 
experimental system or a real metamorphic rock in the 
making - is defined by the values of certain key inten- 
sive variables, including pressure (P), temperature (T) 
and, in multi-component systems, the compositions (X 
values) of one or more phases. For a given equilibrium 
between specific phases, some of these values are auto- 
matically constrained in the phase diagram by the 
equilibrium phase assemblage. The variance of this 
equilibrium is the number of the variables that remain 
free to adopt arbitrary values, which must be deter- 
mined by some other means if the state of the system is 
to be defined completely. 

Phase diagrams in P-T space 



The need to represent phase equilibrium data in visual 
form on a two-dimensional page leads to the use of 
various forms of phase diagram, each having its own 
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merits and limitations. We begin by looking at P-T 
diagrams. 

The two phase diagrams so far considered 
(Figures 1.3a and 2.1) both show the effects of varying 
pressure and temperature on a system consisting of 
only one component (CaC0 3 or AI,SiO-)- Other impor- 
tant examples of such unary systems are discussed in 
Box 2.2. 

P-T diagrams can also be used to show the pres- 
sure-temperature characteristics of multicomponent 
reactions and equilibria. An example is shown in 
Figure 2.2. The univariant boundary in this diagram 
represents not a phase transition between different 
forms of the same compound, but a reaction or equilib- 
rium between a number of different compounds: 




77° C 



NaAlSi,0 6 + SiO z ^ NaAlSi 3 0 8 

jadeite quartz albite /o o\ 

(a pyroxene) (a feldspar) ' ' ' 

For this reason the term reaction boundary (or equilib- 
rium boundary) is used. It marks the P-T threshold 
across which reaction occurs, or the conditions at 
which univariant equilibrium can be established. 

Two components are sufficient to represent all 
possible phases in this system. We can choose them 
in a number of equivalent ways; selecting Na AISi,0, 
and Si0 2 is as good a choice as any. Applying the 
Phase Rule to point X: 



Figure 2.2 P-T diagram showing the experimentally 
determined reaction boundary (solid line) for the reaction 

jadeite + quartz — » albite 

Jadeite (NaAlSi 2 0 6 ) and albite (NaAlSi 3 0 8 ) are both 
aluminosilicates of sodium (Na). 

this system means that a three-phase assemblage is no 
longer invariant, as it was in Figure 2.1. 

At first glance one might expect the albite field 
(point Z) to be divariant like the jadeite + quartz field, 
but here the Phase Rule springs a surprise: 



Point X </i = 2 (2 phases, jadeite + quartz) 

C = 2 (2 components, NaAISi 2 O g and Si0 2 ) 

2+F=2+2 

Therefore F=2 signifying a divariant field. 



Point Z (f> = l (1 phase, albite) 

C=2 (2 components, NaAISi 2 0 6 and Si0 2 ) 

l+F=2+2 

Therefore F= 3 a trivariant equilibrium. 



The jadeite + quartz field is therefore a divariant field 
like that of kyanite in Figure 2.1. At point Y on the 
phase boundary, however, three phases are in equilib- 
rium together: 



Point Y (f> = 3 (3 phases, jadeite + quartz + albite) 

C=2 (2 components, NaAISi 2 0 6 and Si0 2 ) 

3+F=2+2 

Therefore F=1 a univariant equilibrium. 



The three-phase assemblage represents a univariant 
equilibrium: only one variable, P or T, needs to be 
specified to determine completely the physical state of 
the system. The value of the other can be read off the 
reaction boundary. The existence of two components in 



Analysing the albite field in this way, it appears nec- 
essary to specify the values of three variables to 
define the state of the system in this condition. P and 
T account for two of them, but what can the third 
variable be? The answer becomes clear if we ask 
what requirements must be met if, in passing from X 
to Z, we are to generate albite alone. If the mixture of 
jadeite and quartz contains more molecules of Si0 2 
than NaAlSi 2 0 6 , a certain amount of quartz will 
be left over after all the jadeite has been used up. 
The resultant assemblage at Z will therefore be 
albite + quartz. The presence of two phases leads to a 
variance of 2 for this field, as originally expected. 
Conversely, if we react SiO, with an excess of 
NaAlSi 2 0 molecules, the resultant assemblage at Z 
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Box 2.2 Other one-component phase diagrams 



Graphite-diamond 

The phase relations between the two main crystalline 
forms of carbon (Chapter 8) are shown in Figure 2.2.1a. 
Note the very high pressure (in excess of 20xl0 8 Pa) 
required to stabilize diamond at relevant temperatures. 
For this reason, diamond can only form naturally deep 
inside the Earth’s mantle (an equivalent depth scale is 
shown at the right-hand side). Moreover, the minimum 
pressure increases with temperature, so still higher 
pressures are necessary to stabilize diamond in the hot 



interior of the Earth than would be the case at room tem- 
perature. The curve marked ‘geotherm’ shows how tem- 
perature increases with depth beneath ancient continental 
shields (the geological setting where diamond-bearing 
kimberlites are found). From the point at which the geo- 
therm enters the diamond stability field, it is clear that 
diamonds can only be formed at depths greater than about 
120 km (=40 x 10 s Pa). 

Each stability field in Figure 2.2.1a is divariant, and the 
boundary between them is univariant. Because there are 
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Figure 2.2.1 (a) The P-T phase diagram for the two main forms of elemental carbon, graphite and diamond. The curve 
marked ‘geotherm’ shows how temperature increases with depth beneath ancient continental cratons. (b) The P-T 
phase diagram for the system H 2 0 (ice-water-vapour). The inset sketches the corresponding phase diagram for C0 2 . 



would consist of albite+jadeite, again a divariant 
assemblage. The only way to form the one-phase 
assemblage - albite alone - is to combine jadeite and 
quartz in exactly equal molecular proportions, so 
that no quartz or jadeite is left over. In other words, 
to generate just albite in passing from X to Z we must 
control not only P and Tbut also a compositional prop- 
erty of the system - the NaAlSi,0 6 : SiO, ratio. This 



compositional requirement is the unsuspected third 
degree of freedom whose existence the Phase Rule 
has uncovered. 

P-T diagrams provide useful vehicles for portraying 
metamorphic conditions in the crust (Yardley, 1989) 
and for showing the way in which they change with 
time (so-called P-T-t paths) during mountain building 
episodes (Barker, 1998; Best, 2002). P-T diagrams are 
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only two polymorphs of carbon (i.e. </><3), there is no invari- 
ant point (so that F> 0) in this phase diagram. 

The H 2 0 phase diagram 

Figure 2.2.1b shows the phase equilibria between the famil- 
iar forms of pure H 2 0 as a function of pressure and tempera- 
ture. Note that the axes have not been drawn to scale. 
Vapour, liquid and solid coexist at only one point in the 
diagram (0.06xl0 5 Pa and 0.008 °C). This triple point ‘T’ 
lies below atmospheric pressure, which is shown as 
P A =lxl0 5 Pa (dashed line). The curve T-C shows the 
vapour pressure at which liquid water and vapour are in 
mutual equilibrium (the equilibrium or saturation vapour 
pressure) as a function of temperature. Along this curve, the 
vapour is said to be saturated; at pressures below the 
phase boundary T-C, however, the vapour is unsaturated 
and no liquid water can form. The equilibrium vapour pres- 
sure curve T-C rises with temperature to reach P A at 100 °C. 
We can define the boiling point of pure water (T b ) as the 
temperature at which the equilibrium vapour pressure (the 
univariant curve T-C) becomes equal to atmospheric pres- 
sure. The vapour then exerts sufficient pressure to displace 
its atmospheric surroundings and form bubbles in the liquid, 
the everyday phenomenon of boiling. (Note that if the atmos- 
phere pressure is lowered below P A , on a high mountain for 
example, water can boil at temperatures lower than 100 °C.) 

At room temperature (25 °C), the equilibrium vapour pres- 
sure of water lies well below atmospheric pressure P A . 
Atmospheric water vapour can be considered to exert a 
partial vapour pressure in proportion to its concentration in 
the air. If this partial pressure is below the saturation 
vapour pressure, there will be a net evaporation of water 
to vapour (clothes dry, puddles evaporate), whereas if the 
partial pressure of water exceeds the equilibrium vapour 



pressure, water vapour will tend to condense to liquid water 
(dew, mist, rain). The relative humidity expresses the actual 
water vapour pressure in a given body of air as a percent- 
age of the saturation vapour pressure at the temperature 
concerned. 

The univariant curve dividing the liquid and vapour fields 
ends abruptly at invariant point C, known as the critical 
point of water. At this combination of P and T, the structural 
distinction between liquid and gaseous states vanishes. 
The two states merge into a single phase. At higher tem- 
peratures and pressures, H 2 0 exists as a homogeneous 
single phase called a supercritical fluid, in which are com- 
bined the properties of a highly compressed gas and a 
superheated liquid. Some of the hydrothermal fluids respon- 
sible for depositing ore bodies come into this category, and 
the noun ‘fluid’ used alone often has this connotation in 
geology. All liquid/gas systems become supercritical fluids 
under sufficiently extreme conditions. 

The atmospheric pressure (P A ) isobar cuts the ice-water 
phase boundary at exactly 0°C (TJ. Note that this phase 
boundary has a negative slope , 3 a unique feature of the 
ice-water system upon which every ice-skater uncon- 
sciously depends. It expresses the fact that the melting 
point of ice decreases as pressure is increased, so that ice 
close to 0°C can be melted simply by the application of 
pressure, such as the skater’s weight acting on the narrow 
runner of the skate. This behaviour, like many other proper- 
ties of water (Box 4.1), is unique among the common liq- 
uids: melting points for most other materials rise with 
increasing pressure, as the phase diagram for carbon diox- 
ide (inset in Figure 2.2.1b) illustrates. See Exercise 2.2 at 
end of this chapter. 



3 i.e. a negative gradient - see Figure Alb in Appendix A. 



commonly shown inverted, with pressure increasing 
downwards and with a depth scale added to the pres- 
sure axis (cf. Figure 2.5.1). This orientation allows pres- 
sure-dependent phase equilibria to be more easily 
correlated with geophysical profiles of the crust and 
upper mantle. Such diagrams are particularly helpful 
for representing mantle melting processes as a func- 
tion of depth (Box 2.5). 



P-T diagrams 

Reactions in which all the reactants and products are 
crystalline minerals such as those illustrated in 
Figures 1.3a, 2.1 and 2.2 are known as 'solid-solid reac- 
tions'. No vapour is involved, and its presence or 
absence in the experiments is immaterial to the equi- 
librium finally obtained (although it can accelerate 
progress toward equilibrium). Figure 2.3 illustrates 
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Figure 2.3 P H20 - T diagram showing the 'dehydration 
curve' of the mica muscovite (KAl 2 Si 3 AlO 10 (OH) 2 ), showing 
the P H2 0 - T conditions at which it breaks down into the 
assemblage sanidine (KAlSi 3 O g , an alkali feldspar) + 
corundum (A1 2 0 3 )+ vapour (H 2 0). Muscovite and sanidine 
are aluminosilicates of potassium (K). 



another important class of reactions, in which a volat- 
ile constituent plays an essential role. 

It shows a reaction involving water (the dehydration 
of muscovite at high temperatures): 



KAl 3 Si 3 O 10 (OH) 2 ^ KAlSi 3 O g + A1 2 0 3 + H 2 0 

muscovite sanidine corundum vapour 

[arnica) [a feldspar) 



(2.4) 



and Figure 2.1.1a in Box 2.1 shows a similar reaction 
involving carbon dioxide, which is important in the 
metamorphism of siliceous limestones: 



CaC0 3 + Si0 2 ^ CaSiO, +C0 2 

calcite quartz zvollastonite vapour 

(a pyroxene -like (2.5) 

mineral ) 

Because molecules of H,C) and C0 2 are involved in 
these reactions, vapour pressure exerts a strong influ- 
ence on the position of equilibrium. 

The experiments from which these diagrams were pre- 
pared were carried out in the presence of an excess of H 2 0 
or C0 2 respectively (as a separate gas phase), so that at all 
times the experimental charges were saturated with the 



volatile constituent. The consequent existence of a vapour 
in all the experiments means that the vapour pressure 
was equal to the total pressure applied to the specimen. 
The general symbol for vapour pressure is P y ; one can use 
the more specific symbols P H20 or P COz for the vapour 
pressure of water and carbon dioxide respectively. 

Figure 2.3 is therefore a vapour pressure-temperature 
diagram, which can be analysed using the Phase Rule in 
the same way as Figure 2.2. For example, if we choose the 
components carefully in Figure 2.3 (Equation 2.4), we find 
that only three are necessary to constitute all four phases: 

KAlSi 3 0 8 A1 2 0 3 H,0 

since the composition of muscovite can be represented 
as a combination of these three. 



Point F $ = 4 (4 phases = muscovite + sanidine + 

corundum+ water) 

C = 3 (3 components, KAISi 3 0 8 +AI 2 0 3 +H 2 0) 

4+F=3+2 

Therefore F=1 a univariant equilibrium. 



Le Chatelier’s principle 

Behind the empirical facts of mineral stability, as repre- 
sented by the experimentally determined P-T diag- 
rams in Figures 2.1, 2.2 and 2.3, there are some 
important thermodynamic principles which will help 
in the interpretation of phase diagrams. 

The first of these concerns the distribution of phases in 
a phase diagram. Why is kyanite stable at high pressures, 
whereas andalusite can survive only at low pressure 
(Figure 2.1)? What properties of the two minerals dictate 
this behaviour? Flow is it that sillimanite is more stable 
than either of them in the highest temperature range? 

The answers to these questions lie in a simple princi- 
ple enunciated by the French chemist Henri Louis Le 
Chatelier in 1884: when a system at equilibrium experi- 
ences a change in physical conditions, the system will adapt 
in a direction which tends to nullify the change. In the pre- 
sent context, the 'physical conditions' referred to are 
pressure and temperature. 

Consider a system comprising kyanite and andal- 
usite in mutual equilibrium, for example under the 
conditions represented by point D in Ligure 2.1: 

Al 2 SiO s ^ Al 2 SiO s 

andalusite kyanite ( 2 . 6 ) 
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How will the equilibrium respond if we attempt to 
increase the applied pressure (without changing the 
temperature)? Le Chatelier's principle suggests that 
the assemblage will adjust into a more compact form, 
because by taking up less space it can relieve the addi- 
tional pressure applied. The system can accomplish 
this by recrystallizing andalusite (density 3.2 kg dm -3 ) 
into the denser polymorph kyanite (density 3.6 kg dm -3 ). 
By allowing the proportion of kyanite to increase at the 
expense of andalusite, the system can for the time 
being prevent any increase in pressure, hence - in the 
words of Le Chatelier's Principle - 'nullifying the 
change'. Eventually, however, the andalusite becomes 
exhausted and the pressure, no longer constrained by 
univariant equilibrium, is able to rise into the kyanite 
field. From Le Chatelier's principle one can therefore 
show that in any P-T phase diagram the higher-den- 
sity (lower molar volume) phase assemblage will be 
found on the higher-pressure side of a reaction bound- 
ary. Diamond and liquid water are other examples 
(Box 2.2). 




A second consequence of Le Chatelier's rule is that 
the phase assemblage on the high-temperature side of 
an equilibrium boundary is invariably the one having 
the higher enthalpy (Box 4.2). 

The Clapeyron equation 

A second useful application of thermodynamic data to 
phase diagrams is to estimate the gradient (slope) of 
an equilibrium boundary in P-T space. In the symbol- 
ism of calculus (see 'Differentiation' in Appendix A), 

this is written meaning the rate at which P 

dT 

increases for a given increase of T as one follows the 
univariant boundary. The gradient has a sign (positive 
or negative - see Figure Al(b) in Appendix A) and a 
numerical magnitude (indicating whether it is gentle 
or steep). 

Consider the reaction of Equation 2.3 and Figure 2.2: 



The tree-energy change of this reaction is: 

cts ^reactants 

(^jadeite ^quartz ) 



^products ^reactants 



= G., 



(2.7) 



The molar free energy of each phase (which could be 
calculated from published tables of molar enthalpy and 
entropy) varies with pressure and temperature. AG 
therefore varies systematically across a P-T diagram. 
The equilibrium boundary in Figure 2.2 marks the locus 
of P-T coordinates for which AG = 0. It can be shown 
(using quite simple calculus) that the condition for 
remaining on the univariant equilibrium boundary as P 
and T are varied by small amounts dP and dT is that: 



dP AS 
dT ~ AV 



( 2 . 8 ) 



where AS and AV are the entropy and volume changes 
occurring during the reaction in Equation 2.3: 

-Lilbite (^jadeite ^quartz ) 

= 17,!,,^ — ( Vj adeite + V quartz ) 

S and V represent the molar entropy and molar vol- 
ume of each phase, which can be looked up in tables of 
thermodynamic data for minerals (such as Holland 
and Powell, 1998). 

Equation 2.8 is called the Clapeyron equation, after its 
originator, an eminent 19th-century French railway 
engineer. It provides a means of estimating the gradi- 
ent of a reaction boundary in a P-T diagram from eas- 
ily obtainable thermodynamic data. It is also very helpful 
in interpreting many features of phase-equilibrium 
diagrams. 

To predict the slope of the phase boundary in 
Figure 2.2, one proceeds as follows. The relevant molar 
entropy and molar volume data are: 





S 


V 




J K -1 moh 1 


10~ 6 m 3 moh 1 


jadeite (NaAISi 2 O g ) 


133.5 


60.4 


quartz (Si0 2 ) 


41.5 


22.7 


albite (NaAISqOg) 


207.4 


100.1 



Adding the reactants together: 



jadeite + quartz ^ albite 



jadeite + quartz S = 175.0 J K 1 mol 1 V = 83.1 lO^m’mol 1 
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Therefore AS = +32.4 JK -1 mol -1 and AV = +17.0 
10 -6 m 3 mol _1 
Thus: 

dP 32.4 Jmol^K' 1 
df ~ 17.0 10' 6 m 3 mol 1 
= 1.91xl0 6 Jm 3 K 1 
= 1.91xl0 6 Nm 2 K 1 
= 19.1xl0 5 PaK 1 

The units (10 5 PaK _1 =barK _1 ) relate to the gradient of a 
line in P-T space (Figure 2.2). The sign of the gradient 
is positive, consistent with Figure 2.2 (P rises with T), 
and the magnitude (19.1 x 10 5 PaK _1 ) agrees well with 
the value of about 20xl0 5 Pa°C _1 as measured from 
Figure 2.2, noting that kelvins and degrees Celsius are 
units of equal size (Appendix A). 

The positive slope of Figure 2.2 thus reflects the 
observation that both AS and AV for this reaction are 
positive (or both negative, if we write the reaction the 
other way round). A negative slope would signify that 
AS and AV had opposite signs, as is the case for the 
andalusite-sillimanite reaction (Figure 2.1). 

The most striking difference between Figures 1.3a, 
2.1 and 2.2, on the one hand, and Figure 2.3 on the 
other, is that the latter has a curved reaction boundary, 
whereas the others are straight. The reaction boundary 
is curved because the volume change for such a reac- 
tion (and therefore dP/dT) is very pressure-sensitive: 

muscovite sanidine + corundum + vapour 

AV = V + V ... + V , -V 

vapour sarudine corundum muscovite 

At low pressures the volume of the 'vapour' phase 
(actually a supercritical fluid - see Box 2.2) is much 
greater than those of the solid phases, and therefore 
dominates the value of AV. 

AV = V 

vapour 

Because this term is large, the reaction boundary at low 
pressure has only a moderate slope. But the vapour, 
like any gas, is much more compressible than the solid 
phases. At higher pressures, V rapour and AV will get 
progressively less, and the slope of the dehydration 
boundary will get correspondingly steeper. This gen- 
eral shape is a feature of all reactions involving the 
generation of a 'vapour' phase (see also Figure 2.1.1). 
Curved boundaries in P-T diagrams always signify 



the involvement of a highly compressible phase, 
usually a gas (e.g. Figure 2.2.2). 

Phase diagrams in T-X space 



Crystallization in systems with no solid solution 

P-T and P v -T diagrams make no provision for changes 
in the compositions of individual phases during reac- 
tions. Such changes are an important feature of igne- 
ous processes, and make it necessary to introduce 
another type of diagram in which the temperature of 
equilibrium is plotted as a function of phase compos- 
ition ('X'). An example is shown in Figure 2.4, which 
shows the phase relations at atmospheric pressure for 
the binary system CaMgSi 2 0 6 -CaAl 2 Si 2 0 s . Because this 
system is relevant to igneous rocks (it includes simple 
analogues of basalt), the temperature range extends up 
far enough to include melting. 

If the temperature is sufficiently high, it is possible 
to make a homogeneous melt containing the two com- 
ponents CaMgSi 2 O fi and CaAl 2 Si 2 O g in any desired pro- 
portion. These compounds are said to be completely 
miscible in the melt phase. Consequently, in the field 
marked 'melt', only this single phase is stable. In the 
solid state, however, the two components exist as the 
separate phases diopside (ideal composition 
CaMgSi 2 0 6 ) and anorthite (composition CaAlSi,0 8 ): 
there is no stable homogeneous solid of intermediate 
composition. The area below 1274 °C is therefore a 
two-phase field. 

The two areas ABE and ECD are also two-phase 
fields, each representing equilibrium between a melt 
and one of the crystalline phases. To see how, con- 
sider the line xy. This is an isothermal line at a tem- 
perature whose precise value is unimportant (in this 
case it is 1400 °C). We call this a tie-line, because it 
links ('ties') together the compositions of two phases 
which can coexist stably at this temperature, x repre- 
sents the only composition of the melt that can be in 
equilibrium with anorthite (composition y) at 1400 °C; 
it consists of 61% CaAl,Si 2 0 8 and 39% CaMgSi,0 6 . If 
the melt were more CaMgSi,0 6 -rich than this (comp- 
osition Xj for example), it would dissolve anorthite 
crystals and thereby increase the Ca AI 2 Si,O s content 
until equilibrium was reached or until the anorthite 
present had all dissolved. If the liquid had the 
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Figure 2.4 Melting relations in the pseudo-binary system C a M gS i O ,-C a A 1 2 S i 2 O g at atmospheric pressure. The horizontal 
ruling represents two-phase fields: the solid-solid field is ruled more heavily than the solid-melt fields. 'Di' and 'An' refer to 
the phases diopside (composition CaMgSi 2 O s ) and anorthite (composition CaAI,Si 0 O 8 ). This phase diagram is not strictly 
binary because small amounts of aluminium enter the pyroxene phase (see Morse, 1980, pp. 53-7, for details). The circular 
cartoons on the right illustrate what an experimental 'magma' might look like under the microscope at each stage. 



composition x 2 , it would precipitate anorthite, so 
reducing the CaAl 2 Si 2 O g content. 

As the line DE shows, the composition of melt that 
can exist in equilibrium with anorthite (An) depends 
on the temperature. The corresponding line A-E shows 
that the same is true of the melts that can coexist with 
diopside (Di). The curve AED, the locus of the melt 
compositions that can coexist in equilibrium with 
either diopside or anorthite at different temperatures, 
is called the liquidus. All states of the system lying 
above it in Figure 2.4 consist entirely of the melt phase. 
The one point common to both limbs of the liquidus is 
E, which therefore represents the unique combination 
of melt composition and temperature at which all three 
phases are simultaneously at equilibrium. This condi- 
tion is called a eutectic. 

In applying the Phase Rule to Figure 2.4, we must 
recognize that a T-X diagram like this is no more than 
the end-view of more complex phase relations encoun- 
tered in P-T-X space (cf. Figure 2.6a). By considering 
melting relations only at a single - in this case atmos- 
pheric - pressure, we are in fact artificially restricting 
the variance of each equilibrium. Any statements we 



make about variance in this diagram relate only to an 
apparent variance F' where: 

F'= F-l (2.10) 

One may write the Phase Rule in terms of F' as follows: 
</> + F = 0 + (F'+l) = C + 2 
Therefore 

(j) + F’=C + l (2.11) 

This form of the Phase Rule, applicable to isobaric T-X 
(and, incidentally, isothermal P-X) phase diagrams, is 
sometimes known as the Condensed Phase Ride. 



Point E i/) = 3 (3 phases, Di + An + melt) 

C=2 (2 components, CaMgSi 2 0 6 and 

CaAI 2 Si 2 O g ) 

3 + F'=2 + l 

Therefore F'= 0 an isobarically invariant equilibrium 



The term isobarically invariant is jargon that reminds us 
that this equilibrium is invariant only so long as the 
pressure is held constant; if this constraint were to be 



Box 2.3 Tie-lines and the Lever Rule 



In T-X, P-X or P-T-X diagrams, tie-lines link together the 
compositions of two different phases that can coexist in 
equilibrium under specific conditions. Any composition 
lying between the ends of a tie-line must therefore repre- 
sent a physical mixture of the two phases. From the posi- 
tion of that point on the tie-line, one can work out the 
relative proportions of the two phases in the mixture. 

Figure 2.3.1a shows part of a phase diagram in which 
complete solid solution exists between two compounds, A 
and B (cf. Figure 2.5). The tie line c-d depicts equilibrium 
at temperature T 1 between a melt of composition c on the 
liquidus, and a solid solution of composition d on the soli- 
dus; both c and d are expressed in mass % B. Composition 
x lies in the two-phase field between c and d, and must 
signify a physical mixture of these two distinct phases. Let 
C and D represent the mass fractions (i.e. C+D = 1.00) in 
which c and d are mixed to form x. We can then express 
the composition of x as a weighted average of c and d: 

x = Cc + Dd (2.3.1) 

Since C=1 - D, this can be rewritten: 

x=(l - D)c + Dd = c-Dc + Dd 

Therefore x - c = D(d - c) 
x — c 

leading to D = 

d-c 



Substituting D = 1 - C into Equation 2.3.1, we can show in 
a similar fashion that: 



- x-d d-x , , , .. , , 

C = = (thechangeinsigncancelsout) 

c-d d-c 

The mass ratio in which c and d are present in x is there- 
fore given by: 



C _ d - x lx-c 
D d-c i d-c 



(after cancelling denominators) 



Inother words : C(x-c) = D(d-x) (2.3.2) 



This useful equation is known as the Lever Rule, since 
it can also be applied to the ‘lever effect’ of the old- 
fashioned beam-balance (Figure 2.3.1b), in which the 
weight of a body C is inversely proportional to the dis- 
tance from the fulcrum (c-x) at which it balances an 
opposing weight D: 



weight of C _ x-d 
weight of D c-x 



Qualitatively, the closer the composition of a mixture 
plots (in composition space) to one of its constituents, 
the greater the percentage of that constituent in the 
mixture. 




Figure 2.3.1 (a) Part of a phase diagram similar to Figure 2.5 to illustrate the Lever Rule, (b) The analogous geometry 
of the beam-balance. 
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relaxed, we would find that the eutectic was just one 
point of a univariant curve in T-X-P space. 

A eutectic is always an invariant point in any phase 
diagram. On the other hand, the one-phase 'melt' field 
has two degrees of freedom: 



Point x 1 (/> = 1 (1 phase, melt) 

C= 2 (2 components, CaMgSi 2 0 6 and 

CaAI 2 Si 2 0 8 ) 

1 + F'=2 + 1 

Therefore F'= 2 an isobarically divariant equilibrium 



T and X have to be quantified to define fully the state 
of the system in this condition. 

It would be natural to expect the fields ECD and 
ABE to be divariant as well, but the Phase Rule indi- 
cates otherwise. Consider the composition represented 
by point x r Neither the melt (at this temperature) nor 
anorthite can have this composition. The 'composition' 
x 2 only has meaning at 1400 °C when interpreted as the 
composition of a physical mixture of melt x and anor- 
thite y. (The proportions of the two phases in this mix- 
ture can be worked out as explained in Box 2.3.) ECD 
and ABE are therefore two-phase fields. If </>=2 and 
C = 2, we cannot escape the conclusion that F' = l. In 
other words, specifying temperature is sufficient to 
define the composition of all phases in equilibrium, or 
vice versa given that the pressure is already defined. 
Thus reaction-boundary lines (cf. Figures 1.3a and 2.1) 
are not the only manifestation of univariant equilib- 
rium in phase diagrams; areas can also be univariant. 
Such fields arise in T-X diagrams whenever two coex- 
isting phases have different compositions. One can 
imagine them consisting of an infinity of horizontal tie- 
lines, as the horizontal ruling in Figure 2.4 is intended 
to symbolize. 

T-X diagrams are important in igneous petrology 
because they allow one to follow the evolution of melt 
composition with advancing crystallization in experi- 
mental and natural magmatic systems (at constant 
pressure). Imagine a melt m cooling from some temp- 
erature above the liquidus, say 1450 °C. The phase 
assemblage at this point consists of melt alone, as illus- 
trated by the top 'microscope view' on the right. At 
first there will be no change other than a fall in tem- 
perature: we can imagine point m falling vertically 
through the 'melt' field. Arrival at the liquidus (point 
x) signals the first appearance of solid anorthite, which 



here begins to crystallize in equilibrium with the melt. 
Extraction of anorthite depletes the melt a little in the 
CaAl 2 Si 2 O g component, causing a shift in composition 
to the left in Figure 2.4. But the maintenance of univari- 
ant equilibrium demands that the melt composition 
should change in conjunction with falling tempera- 
ture. As crystallization advances, therefore, the melt 
composition migrates steadily down the liquidus 
curve towards E, continuously crystallizing anorthite 
and changing composition. The phase assemblage here 
is shown in the middle roundel on the right, illustrat- 
ing what might be seen in a quenched sample under 
the microscope. 

On reaching the eutectic, the melt begins to crystal- 
lize diopside in addition to anorthite, as shown in the 
lowermost roundel. At this juncture, the melt composi- 
tion becomes fixed, because anorthite and diopside 
crystallize in the same proportion as the CaAI,Si,O g : 
CaMgSi 2 0 6 ratio of the melt. The temperature also 
remains constant, because the eutectic is an invariant 
equilibrium (at least within the isobaric framework 
being considered): as long as three phases remain in 
equilibrium, neither melt composition nor tempera- 
ture can change. Continued cooling in this context 
merely means the loss of heat (the latent heat of crys- 
tallization of diopside and anorthite) from the system 
at constant temperature and the formation of crystals 
at the expense of melt. Eventually the melt becomes 
exhausted, and invariant (Di + An + melt) equilibrium 
gives way to univariant (Di + An) equilibrium, allow- 
ing the temperature to resume its downward progress. 
The total solid assemblage will obviously now consist 
of 61% anorthite and 39% diopside (c in Figure 2.4). 

The eutectic therefore represents the lowest projec- 
tion of the melt field, the composition and temperature 
of the last melt to survive during the cooling of the sys- 
tem. Progressive crystallization of any melt in this sys- 
tem (with the special exceptions of pure CaMgSi 2 0 6 
and pure Ca AI ,Si,O s ) will lead its composition ulti- 
mately to the eutectic. This illustrates an important 
general principle of petrology: that the evolving com- 
positions of crystallizing magmas of all types tend to 
converge upon one or two 'residual magma' composi- 
tions (a natural example being granite) at which the 
liquidus reaches its lowest temperature. 

The eutectic also indicates the composition of the 
first melt to appear upon heating any mixture of diop- 
side and anorthite (Box 2.4). 
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Box 2.4 Partial melting I: melting in the laboratory 



Phase diagrams provide many insights into the important 
process of rock melting. Consider first the progressive 
melting of a mixture of diopside and anorthite, for example 
the mixture of composition c in Figure 2.4. 

The mixture will simply heat up until the temperature 
reaches 1274°C. At this point a melt of the eutectic com- 
position E begins to form; this is the only melt composi- 
tion that can be in equilibrium with diopside and anorthite, 
both of which are at this stage present in the solid mix- 
ture. Continued heating brings an increase in the propor- 
tion of melt at constant temperature, with no change of 
melt composition (invariant equilibrium) until the diopside 
disappears, having been entirely incorporated in the melt. 
(Anorthite has been dissolving too, but has not yet been 
used up.) At this stage, from the Lever Rule (Box 2.3), the 
ratio of melt:anorthite is XC:EX. Univariant equilibrium 
(An + melt) now obtains, and with increasing temperature 
the proportion of melt continues to rise, its composition 
proceeding up the liquidus curve as more and more anor- 
thite dissolves in it. At x, the melt has the same composi- 
tion as the solid starting mixture, and here the last 
remaining crystals of anorthite disappear. The system now 
enters the divariant melt field, where the temperature can 
continue rising without further change of state. 

Similar principles govern the melting of solid-solution 
minerals. The olivine system shown in Figure 2.4.1 (analo- 
gous to Figure 2.5) provides an example relevant to basalt 
production by partial melting in the upper mantle (which 
consists largely of olivine). When olivine is heated up to 
the solidus (for example, point c 1 ), a small proportion of 
melt m 1 appears. The melt is much less magnesium-rich 



than the olivine from which it is produced, a point of great 
petrological significance. Continued heating will cause the 
temperature to rise, the proportion of melt to increase, its 
composition to migrate up the liquidus curve towards m 2 , 
and that of the remaining olivine crystals to migrate up the 
solidus towards c 2 . The system would become completely 
molten at just over 1800 °C (m 3 ). Thus a gap exists 
between the temperature at which olivine begins to melt 
and that at which it becomes completely liquid (as in 
Figure 2.5). This gap, called the melting interval, is a fea- 
ture of all minerals that exhibit solid solution. The every- 
day notion of a ‘melting point’ applies only to pure 
end-members, where the liquidus and solidus converge. 

The complete melting of rocks like this only occurs in 
very unusual circumstances, like meteorite impacts. 
Generally, magmas are produced by a process of partial 
melting, in which temperatures are sufficient to melt a 
fraction of the source material but not all of it. Both in the 
olivine phase diagram and in actual rocks, partial melting 
will generate a melt less magnesium-rich (m 2 ) than the 
source material (c 1 ), leaving behind a refractory solid resi- 
due which contains more magnesium (c 2 ) than the source 
material prior to melting. The composition of both prod- 
ucts - melt and residual solid - depends upon the degree 
(percentage) of melting, and therefore on the temperature 
attained. 

We must be careful not to assume that the tempera- 
tures shown in the olivine diagram are necessarily charac- 
teristic of the upper mantle. We have seen that although 
pure anorthite (Figure 2.4) remains solid up to 1553 °C, a 
mixture of anorthite with diopside begins to melt below 



Eutectics are common features in systems of this 
kind. They represent the general observation that mix- 
tures of minerals (in other words, rocks) begin to melt at 
lower temperatures than any of the pure constituents 
(minerals) would on their own, just as a mixture of ice 
and salt conveniently melts at lower temperature than 
ice alone. This principle is widely used in industry 
when a flux is added to enable a substance to melt at a 
lower temperature than it would in the pure state (e.g. 
in soldering). 



Crystallization in systems with solid solution 

Diopside and anorthite belong to different mineral 
groups having different crystal structures, and the ten- 
dency for either to incorporate the constituents of the 
other into its crystal structure is negligible. But within 
many mineral groups it is common to find that crystal 
composition can vary continuously between one end- 
member composition and another. One can visualize 
one solid end-member 'dissolving' in the other to form 
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M 92 Si °4 Percent Mg 2 Si0 4 by mass Fe 2 Si0 4 

Fo component Fa component 

Figure 2.4.1 Melting relations in the olivine series at atmospheric pressure. This system has been determined 
experimentally only up to 1500 °C, and for pure forsterite. Interpolated boundaries are shown dashed. 



1300 °C. The same is true of mantle rocks, which consist 
not of olivine alone (although it is the dominant constit- 
uent), but contain pyroxene and either garnet or spinel as 
well. Any such mixture will begin melting at temperatures 
below those at which the individual minerals would start to 
melt. Moreover, every mineral contributes to the compos- 
ition of the melt right from the start, just as both anorthite 



and diopside contribute to the eutectic melt in Figure 2.4. 
It is a common misconception to believe that minerals in a 
rock will melt one after the other. A partial melt should be 
seen as a solution in which all the solid phases of the 
source rock are partly soluble. 

The influence of pressure (depth) on melting is consid- 
ered in Box 2.5. 



a homogeneous crystal of intermediate composition. 
This phenomenon of miscibility in the solid state is 
referred to as solid solution. 

Figure 2.5 shows the crystallization behaviour of a 
familiar example of such solid solution series, plagio- 
clase feldspar (a solid solution between albite, 
NaAlSi 3 O g , and anorthite, CaAl 2 Si 2 O g ). Only one solid 
phase appears in the diagram. 

Prominent in the diagram is a leaf-shaped feature, 
bounded by two curves depicting coexisting crystal 
and melt compositions as a function of temperature. 



Thus line ab in Figure 2.5 is a tie-line linking melt com- 
position a to the plagioclase crystal composition b with 
which it is in equilibrium at that temperature. The 
curve through a, above which (in the 'melt' field) the 
system is entirely molten, is the liquidus; the curve 
through b, below which the system consists entirely of 
crystalline plagioclase, is called the solidus. 

By applying the Condensed Phase Rule to Figure 2.5, 
we find that both the 'melt' and 'plagioclase' fields are 
divariant: </> = 1, C = 2 (NaAlSi 3 O g + CaAl 2 Si,O g ), there- 
fore F' = 2. The 'melt + plagioclase' field is another 



Box 2.5 Partial melting II: melting in the Earth’s mantle 



Figure 2.5.1 shows in simple terms how the solidus 
temperature of peridotite - the temperature at which it 
begins to melt - varies with depth in the mantle. 4 This is 
in essence a P-T diagram, drawn ‘upside-down’ in terms of 
pressure in order to depict temperature as a function of 
depth below the surface. The band between the solidus 
and liquidus lines, the ‘melting interval’, depicts the range 
of conditions under which partial melting (Box 2.4) of man- 
tle peridotite can occur. These are the conditions neces- 
sary to produce basaltic magma. 

The curve marked ‘Oceanic geotherm’ indicates how the 
ambient temperature in the upper mantle is believed to vary 
with depth beneath a typical sector of mid-ocean ridge. Note 
that the geotherm fails to reach the solidus at any depth. 
Why then should melting occur at all in the upper mantle? 

To see why, we must recognize that the mantle is not a 
static body. Because the interior of the Earth is hot, the 
solid mantle undergoes continuous (although very slow) 
convective motion, with mushroom-like ‘plumes’ of buoy- 
ant hotter material ascending from below (e.g. beneath 
Hawaii), and dense colder material sinking down (e.g. cold 
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oceanic lithosphere at subduction zones). In a converting 
mantle, melting may arise purely as a consequence of 
upward motion. In Figure 2.5.1a, the solid peridotite at 
point X, for instance, can penetrate the solidus and begin 
to melt simply by migrating upward to lower pressures along 
the path X-Y. This process, known as decompression 
melting, is the primary cause of magma generation at 
mid-ocean ridges: plate forces continuously pull the lith- 
ospheric plates apart and thereby permit passive upwelling 
(Figure 2.5.1b) and melting of the underlying astheno- 
sphere. No rise in temperature is required; indeed, the 
ascending material cools slightly as a result of the work it 
has to do in expanding (Box 1.1). 

Mantle plumes, on the other hand, are sites of buoyant 
upwelling of deeper material that may be 150-300 °C hot- 
ter than the surrounding upper mantle. Melting in a plume 
is the combined effect of an elevated geotherm and 
decompression resulting from upwelling. As many plumes 
are located in intra-plate settings, the presence of thick 
cool lithosphere confines melting to deeper levels than 
beneath a mid-ocean ridge. 
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Figure 2.5.1 (a) Variation with pressure and depth 
of the solidus temperature of dry mantle peridotite. 
The heavy curve labelled ‘oceanic geotherm’ shows 
how ambient temperature beneath an oceanic ridge 
is believed to vary with depth. X-Y illustrates a 
passive upwelling path resulting from lithosphere 
extension, (b) Cartoon illustrating where 
decompression melting occurs beneath an exten- 
sional mid-ocean ridge and in an intraplate mantle 
plume head. 
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4 In the absence of water. 
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NaAISi 8 0 8 CaAI 2 Si 2 0 8 

(Ab component) (An component) 

Percent CaAI 2 Si 2 0 8 
by mass 

Figure 2.5 T-X diagram showing melting relations in 
plagioclase feldspars at atmospheric pressure. Horizontal 
ruling represents a two-phase field. The subscript 'ss' 
denotes that plagioclase is a solid solution (see text). 

example of a univariant field (cf. Figure 2.4): t/> = 2,C = 1, 
therefore F' = 1. If equilibrium exists between melt and 
plagioclase, specifying T automatically defines the 
compositions of both phases. Conversely, knowing 
either phase composition defines the other composi- 
tion and the temperature of equilibrium unambigu- 
ously. The line ab is one of an infinite series of such 
tie-lines traversing the two-phase field, as the horizon- 
tal ruling symbolizes. 

Crystallization in this system leads to a series of con- 
tinuously changing melt and solid compositions. Melt 
m, for example, will cool until it encounters the liquidus 
curve at a, where plagioclase b will begin to crystallize. 
Because b is more CaAl,Si,0 8 -rich than a, its extraction 
will deplete the melt in CaAl 2 Si,0 6 and thereby enrich it 
in NaAlSi 3 O g ; with continued cooling and crystalliz- 
ation the melt composition will migrate down the liqui- 
dus curve. The changing melt composition causes a 
corresponding evolution in the equilibrium composi- 
tion of the plagioclase crystals. Not only will newly 
crystallized plagioclase be more albite-rich than b, but 
there will be a tendency for early-formed crystals, by 



continuous exchange of Na, Ca, A1 and Si, to re-equili- 
brate with later, more albite-rich fractions of the melt. 
Thus to maintain complete equilibrium as the melt 
evolves to point a 1 on the liquidus, all the crystals that 
have so far accumulated must adjust their composition 
to b 1 on the solidus. Such adjustment requires solid-state 
diffusion to and from the centre of each crystal, and is a 
slow process. Crystal growth during natural magma 
crystallization commonly proceeds too quickly to allow 
complete continuous equilibrium between crystals and 
magma: the surface layers readjust to changing melt 
composition, but the crystal interiors fall behind. The 
result is a compositional gradient between relatively 
anorthite-rich cores and relatively albite-rich margins of 
the plagioclase crystals, a phenomenon known among 
petrologists as zoning (Box 3.1). Zoning is seen in other 
mineral groups as well, notably pyroxene. 

If crystallization proceeds slowly enough to permit 
continuous re-equilibration (an ideal situation known 
as equilibrium crystallization), the final melt has the 
composition a 2 and the end-product is a mass of crys- 
tals all having the composition b , the same as the 
initial melt composition m. Imperfect re-equilibration 
between crystals and the evolving melt, however, ties 
up a disproportionate amount of the anorthite compo- 
nent in the cores of early formed crystals owing to 
overgrowth or burial, and the melt can then evolve to 
compositions beyond a n before it runs out, creating late 
zones of plagioclase more albitic than the original melt. 
This process, in which isolation of early-formed solids 
allows later melts to develop to more extreme compo- 
sitions, is called fractional crystallization. Crystallization 
of natural melts in crustal magma chambers approxi- 
mates closely to fractional crystallization, and this pro- 
cess contributes a lot to the chemical diversity of 
igneous rocks and magmas. 

The solvus and exsolution 

The final T-X section to be examined (Figure 2.6b) 
shows the phase relations of the alkali feldspars in the 
presence of water vapour (at a pressure of 2 x10 s 
Pa = 2 kbars). This diagram can be visualized as a cross- 
section of P H , 0 - T-X space, coinciding with the plane 
in which P HO = 2x 10 8 Pa (Figure 2.6a). One can speak 
of the diagram as an isobaric T-X section of phase 
relations in P H2Q - T-X space. 



Box 2.6 Reaction points and incongruent melting 



Every geology student knows that olivine and quartz are 
incompatible, and do not coexist stably in nature. (In fact, 
this is true only of magnesium-rich forsteritic olivines. 
Fayalite - Fe 2 Si0 4 - is quite a common mineral in granites 
and quartz syenites.) How is this incompatibility expressed 
in a phase diagram? 

The relevant part of the system Mg 2 Si0 4 -Si0 2 (omitting 
complications at the Si0 2 -rich end) is shown for atmos- 
pheric pressure in Figure 2.6.1. In many respects it is sim- 
ilar to Figure 2.4. The difference is that, between Mg 2 Si0 4 
and Si0 2 along the composition axis, lies the composition 
of the pyroxene enstatite, Mg 2 Si 2 0 6 . Consider the crystal- 
lization of melt composition m r On reaching the liquidus it 
will begin to crystallize olivine, whereupon further cooling 
and crystallization will lead the melt composition down the 
liquidus curve. On reaching R, the melt composition has 
become too Si0 2 -rich (more so than enstatite) to coexist 
stably with olivine, which therefore reacts with the Si0 2 in 
the melt to form crystals of enstatite: 

Mg 2 Si0 4 + Si0 2 -> Mg 2 Si0 6 

olivine melt pyroxene ( 2 . 6 . 1 ) 

(This symbolism does not mean that the melt consists of 
Si0 2 alone. Other components are present, but this reac- 
tion involves only the Si0 2 component.) 

At R, the three phases are at equilibrium. Using the 
Condensed Phase Rule, it is clear that R is an invariant 



point like E. It is called a reaction point. Temperature and 
melt composition remain constant as the reaction pro- 
ceeds (from left to right in reaction 2.6.1), until one or 
other phase is exhausted. In this case (beginning with m 4 ) 
the melt is used up first, and the final result is a mixture 
of olivine and enstatite: the melt never makes it to the 
eutectic. If, on the other hand, the initial melt had the 
composition m 2 , more siliceous than enstatite, the reac- 
tion at R would transform all of the olivine into enstatite, 
with some melt left over. The disappearance of olivine 
releases the system from invariant equilibrium R, and the 
melt can proceed down the remaining liquidus curve, crys- 
tallizing enstatite directly until the eutectic is reached. The 
final result is a mixture of enstatite and silica (the high- 
temperature polymorph cristobalite). The proportions in 
the final mixture can be worked out by applying the Lever 
Rule to m 2 . 

During melting, this reaction relationship manifests 
itself as a phenomenon called incongruent melting. 
Pure enstatite, when heated, does not melt like olivine 
or anorthite but decomposes at 1557 °C to form olivine 
(less Si0 2 -rich) and melt (more Si0 2 -rich than itself), 
i.e. the reaction 2.6.1 run in reverse. The system is 
held in invariant three-phase equilibrium until the 
enstatite has been exhausted, then continues melting 
by progressive incorporation of olivine into the melt 
(cf. Figure 2.4). 
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Figure 2.6.1 The Mg-rich part of the system 
Mg 2 Si0 4 -Si0 2 showing the reaction point R 
between Mg 2 Si0 4 (forsterite) and Si0 2 -rich 
melts (reaction 2.6.1). 

' J 



EQUILIBRIUM IN GEOLOGICAL SYSTEMS 



35 



(a) 




(b) 





(Or) by mass (Ab) 



Figure 2.6 Melting and subsolidus phase relations in the alkali feldspars (the system KAlSi 3 0 8 -NaAlSi 3 0 g ). The subscript 'ss' 
denotes solid solution, (a) Perspective sketch of P Hi0 - T - X space, showing the isobaric section at 2 x 10 8 Pa illustrated in (b). 
(b) Alkali feldspar phase relations at P HO = 2 x 10 s Pa. Horizontal ruling represents two-phase fields. 



The diagram shows the liquidus and solidus of the 
alkali feldspar series, which differ from those in 
Figure 2.5 only in that they fall to a minimum melting 
point in the middle of the series, rather than at one 
end. As a result we get two leaf-shaped fields instead 
of one. But interest is mainly in what happens in the 
subsolidus region. The 'homogeneous feldspar ss ' field 
immediately below the solidus means that here the 
end-members are completely miscible in the solid 
state: they form a complete solid solution in which 
any composition can exist as a single homogeneous 
phase. But at lower temperatures things get more 
complicated. 

Beneath a boundary called the solvus a 'miscibility 
gap' appears. At these temperatures (for example 
600 °C) the albite crystal structure is less tolerant of the 
KAlSijOg component (partly because of the large size 
of the potassium atom), and at a KAlSi 3 0 8 content of 
about 20% (/, = 80% NaAlSi 3 O g ) becomes saturated 
with it. Any KAlSi 3 O g present beyond this limit is 
forced to exist as a separate KAlSi 3 O g -rich feldspar 
phase, whose composition can be found by extending 
a tie-line to the left-hand limb of the solvus, cutting it 
at /, (about 65% KAlSi 3 O g or 35% NaAlSi 3 0 8 ). This 
potassium feldspar is itself saturated with NaAlSi 3 O g . 



A homogeneous alkali-feldspar solid solution such 
as f h ceases to be stable as it cools through the solvus. At 
point /, for example, it is well inside the two-phase 
region. Such a point represents, at equilibrium, a mix- 
ture of two phases. The initially homogeneous feldspar 
therefore breaks down, or exsolves, into two separate 
phases f 1 and f r But solid-state diffusion is too slow to 
allow a cooling feldspar crystal to sort itself out into 
two separate crystals. The usual result of exsolution is 
a series of thin, lamellar domains of one phase enclosed 
within a host crystal of the other. The lever rule (Box 2.3) 
tells us that in the present example _/j will be more 
abundant than f, and the cooling of crystal/ h will there- 
fore produce a host crystal of composition^ containing 
exsolution lamellae of phase / . Such structures are char- 
acteristic of alkali feldspars, where they are known as 
perthites. Figure 2.7 shows a crystal of perthite viewed 
in a polarizing microscope configured to highlight this 
texture. The dark streaks are albite lamellae; the lighter 
host is orthoclase (divided into upper and lower por- 
tions that differ in shade owing to the twinning of the 
crystal). Exsolution textures analogous to perthite 
(although not given this name) are developed in some 
pyroxenes (Plate 3), owing to a similar miscibility gap 
between diopside (CaMgSi 2 0 6 ) and enstatite (Mg 2 Si 2 0 6 ). 
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Figure 2.7 Part of a twinned crystal of perthite (orthoclase 
host enclosing darker exsolution lamellae of albite) seen 
under a polarizing microscope with crossed polars; width 
of field~ 0.7 mm. 

Application of the condensed phase rule (because 
P lli() is constant) to Figure 2.6b indicates that the sub- 
solvus region is a univariant area. In favourable cir- 
cumstances, coexisting alkali feldspar (or pyroxene) 
compositions can be used to estimate the temperature 
of equilibration, a fact that has practical application in 
geothermometry. 

Ternary phase diagrams 



Merely to represent the range of compositions possible 
in a three-component system requires the use of the 
two dimensions of a piece of paper (Box 2.7). If we 
wish to represent comprehensively the phase relations 
of such a system over a range of temperature (in a form 
analogous to Figure 2.4), we would need to build a 
three-dimensional model. However, as a means of 
disseminating phase equilibrium data, such models 
would be less than convenient, so petrologists have 
devised various ways of condensing their content into 
two-dimensional form that can be printed on a page. 
Examples are shown in Figures 2.8 and 2.9, both of 
which relate to the crystallization of simple silicate 
melts, that is, to laboratory analogues of magmatic 
systems. 

The base of any ternary phase diagram is an equilat- 
eral triangle, in which any composition in a ternary 
system can be plotted (Box 2.7). The temperature axis 
of the hypothetical model is constructed perpendicular 
to the plane of this triangle (see inset to Figure 2.8). 



Ternary phase diagram with no solid solution 

Figure 2.8 shows an example of a ternary phase dia- 
gram involving end-members that exhibit no mutual 
solid solution. The liquidus forms a 3-dimensional 
surface in this diagram, analogous to the liquidus curve 
in Figure 2.4. The surface takes the form of several 
curving 'hillsides' that meet along 'thermal valleys'. 
They can be portrayed either by a three-dimensional 
model constructed on a triangular base - see the inset 
in Figure 2.8 - or by plotting the model in the form of a 
two-dimensional 'map' with temperature 'contours' 
(isotherms) as shown in the main diagram. Note to 
begin with that the left-hand face of the model is equiva- 
lent to the Di- An binary phase diagram that has already 
been discussed in Figure 2.4, albeit shown here in a 
back-to-front sense. In the main diagram the liquidus 
surface consists of four sloping fields separated by 
gently curved boundaries. Each field represents a 
domain of melt composition within which one partic- 
ular mineral crystallizes first (analogous to each curved 
line in Figure 2.4). The largest field is that in which for- 
sterite is this 'liquidus phase'; the fields where diopside 
and anorthite crystallize first are significantly smaller. 

Consider the course of crystallization followed by 
the initial melt composition labelled x. If at the outset it 
has a temperature above the liquidus (say 1500 °C), the 
first stage of its cooling path simply involves cooling 
to the liquidus temperature (about 1430 °C). There, 
because this composition lies within the forsterite field, 
olivine will begin to crystallize (see cartoon (a)). As no 
iron is present in this system, the olivine that forms 
will be pure forsterite (Mg,Si0 4 ). Separating out forst- 
erite in crystalline form depletes the remaining melt in 
the Mg 2 Si0 4 component, a fact that we can represent on 
this diagram by drawing a line (dashed) from the 
Mg,Si0 4 apex to the co-ordinates of x, then extending it 
beyond x (solid arrow). As cooling continues, the melt 
composition will migrate along this arrow, directly 
away from the forsterite apex. Eventually it reaches the 
boundary between the forsterite and diopside fields in 
Figure 2.8. At this boundary, called a cotectic, diopside 
begins to crystallize alongside forsterite. The crystal- 
line extract that is removed from the melt from this 
point on will be a mixture of Di and Fo as in cartoon (b), 
the composition of which must lie somewhere on the 
bottom edge of the diagram (the domain of An-free Di- 
Fo mixtures), and extracting it will therefore drive the 
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Box 2.7 How a ternary diagram works 

Three variables can be plotted in a two-dimensional graph 
if they add up to 100%. Any composition in a three-compo- 
nent system can therefore be represented in two dimen- 
sions, usually in the form of a ternary diagram plotted on 
special equilateral-triangle graph paper (Figure 2.7.1a). 
The user labels each apex with one of the components, as 
shown in Figure 2.7.1b. 

Each apex represents 100% of the component with 
which it is labelled (the Di corner, for example, represents 
a composition consisting of 100% diopside). The side 
opposite represents compositions that are devoid of that 
component (in this case, a range of Ab-An mixtures). 
Lines parallel to this edge are contours representing differ- 
ent Di percentages from 0 (on the Ab-An edge) to 100% 
(at the Di apex). 



To plot a composition such as Di 72%, Ab 19%, An 9% 
(note that the three co-ordinates must add up to 100%) 5 , 
rule a line horizontally across the diagram at the position 
equivalent to 72% Di. Rule another line, parallel to the 
Di-An edge, at the position corresponding to 19% Ab (be 
sure to count this from the Di-An edge where Ab = 0%). 
The intersection with the first line marks the composition 
being plotted. Note that only two readings need to be plot- 
ted; the third - being the difference from 100 - is not an 
independent variable, but it is useful to read it off the 
diagram (An = 9%) to check that the point has been plotted 
accurately. 

As with any plot showing mineral proportions, it is 
essential to indicate whether the numbers relate to mass 
% or molar proportions. 




Figure 2.7.1 (a) Example of ternary graph paper with 100 mm sides (may be photocopied as needed), (b) Flow to plot a 
three-component mixture in a ternary diagram. 



5 



The outcome of plotting three numbers that add up to less than 100% is a smaller triangle surrounding the target point, as illustrated 
in Gill (2010), Fig. Bl, p. 363. In such cases, each value should be multiplied by 100/Tot (where Tot represents the initial total), 
thereby scaling them up to values that do sum to 100%. 

) 
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Figure 2.8 Phase diagram for the ternary system CaAl 2 Sy0 8 -CaMgSi,0 6 -Mg 2 Si0 4 at atmospheric pressure, after Osborne 
and Tait (1952). The liquidus surface is represented as isotherms (^temperature 'contours' graduated in °C). V- ticks mark 
10% graduations in mass % of each component. V and 'y' are illustrative initial melt compositions discussed in the text. 
Circular cartoons represent the 'phenocryst' mineralogy at key stages along each crystallization path, as seen under a 
microscope. E represents a eutectic point (see text) and R represents a reaction point (cf. Box 2.6). The inset sketch from Gill 
(2010) shows how the liquidus would appear as a 3D model (with temperature forming the vertical axis). 



melt composition toward the An apex. In fact, because 
the Fo-Di boundary is a 'thermal valley', the melt com- 
position migrates with falling temperature along the 
cotectic toward the point E - changing melt composi- 
tion in any other direction would require a rise in its 
temperature. Continued removal of Di + Fo eventually 
drives the melt composition to the point E where the 
three fields meet. Flere anorthite begins to crystallize 
alongside the forsterite and diopside (cartoon (c)). 

Melt y will follow a different path as it crystallizes but 
will arrive at the same final melt composition. Initial 
crystallization of forsterite drives melt y to intersect the 
An-Fo cotectic, at which point anorthite will begin to 
crystallize alongside forsterite (see cartoon (d)). Diopside 
will only appear when the melt has reached point E 
(cartoon (c)). This point represents (a) the composition 
at which the liquidus surface reaches its lowest tem- 
perature, and (b) the composition toward which all melts 



will converge as they crystallize, even melts whose 
initial compositions lie in the anorthite or diopside 
fields. It is called the ternary eutectic of this system. 

Applying the Phase Rule to the composition y in 
Figure 2.8 at a temperature above the liquidus yields: 



Point y 

^"liquidus) 


0 = 1 


melt alone 


C=3 

1 + F'= 3 + 1 


CaAI 2 Si 2 0 8 + Mg 2 Si0 4 + CaMgSi 2 0 6 




Therefore 


F'= 3 


a trivariant equilibrium 



The three degrees of freedom calculated here are 
(a) the two compositional coordinates required to 
define the position of y in the ternary diagram, plus (b) 
temperature, which (like point m in Figure 2.4) is 
unconstrained as long as the temperature remains 
above the liquidus. 
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When the temperature falls to meet the liquidus 
(about 1430 °C at point y), the variance of the equilib- 
rium changes as olivine begins to crystallize: 



Point y 


^ = 2 


melt + forsterite 


(T=T -h ) 

' liquidus 7 


C=3 

2+F=3+l 


CaAI 2 Si 2 0 8 +Mg 2 Si0 4 +CaMgSi 2 0 6 


Therefore 


F'= 2 


a divariant equilibrium 



Melt composition y can coexist with olivine only at the 
liquidus temperature, so temperature (now dictated by 
that of the liquidus surface) ceases to be an independ- 
ent variable and the variance decreases to 2. 

When the melt evolves to point z, plagioclase begins 



to crystallize too: 


Point z 


c/> = 3 


melt + forsterite + anorthite 




C = 3 


CaAI 2 Si,0 8 + Mg,Si0 4 + CaMgSi,0 K 




3 + F'=3 + l 




Therefore 


F'= 1 


a univariant equilibrium 



For this equilibrium to be maintained, the two compo- 
sitional parameters can only vary in an interdependent 
way that confines z to the cotectic line, making this a 
univariant equilibrium. 

When crystallization of forsterite and anorthite has 
driven the residual melt composition to point E: 



Point E (/> = 4 melt+forsterite + anorthite + diopside 

C=3 CaAI 2 Si 2 0 8 + Mg 2 Si0 4 + CaMgSi 2 0 6 

4 • r-3 • 1 

Therefore F'= 0 an isobarically invariant equilibrium 



Flere, with the melt in equilibrium with forsterite, 
anorthite and diopside, we have reached an invariant 
situation, the ternary analogue of the binary eutectic in 
Figure 2.4. Melt composition and temperature now 
remain fixed, and only the proportions of the various 
phases can vary: as heat is lost at a constant tempera- 
ture of 1270 °C, melt crystallizes into forsterite, anor- 
thite and diopside until no melt remains. This removes 
one of the four phases (melt) from consideration, leav- 
ing the three solid phases to continue cooling in a 
solid-state, univariant equilibrium. 

A second invariant point R exists in this diagram, 
a reaction point similar to that discussed in Box 2.6. 
In the present context this complication can be 
ignored. 



Ternary phase diagram with solid solution 

Figure 2.9 shows another example of a ternary phase 
diagram relevant to magma crystallization. It incorp- 
orates along its edges two binary phase diagrams 
we have already discussed. Figures 2.4 and 2.5. The 
involvement of the plagioclase series (Figure 2.5) 
introduces solid solution, giving this diagram a dif- 
ferent appearance to Figure 2.8. The gross features of 
the system can again be appreciated from a perspec- 
tive sketch (Figure 2.9a) showing the topography of 
the liquidus surfaces as in Figure 2.8. In this diagram, 
however, these liquidus surfaces meet in a V-shaped 
low-temperature trough running out of the binary 
eutectic in the system Ca MgS i ,0,-Ca A I ,S i ,O g . The 
phase relations in this and the companion binary 
systems (CaMgSi n 0 6 -NaAlSi 3 0 g and CaAl 2 Si 2 O g - 
NaAlSLO.) can be indicated on the vertical faces of 
the 'model'. 

The main diagram (Figure 2.9b) is invaluable for 
examining the evolution of melt composition during 
crystallization, and considering the parallel magmatic 
evolution in real igneous rocks. The V-shaped valley 
divides the diagram into two fields, each labelled with 
the name of the solid phase that crystallizes first from 
melts whose compositions lie within that field. For 
example, a melt of composition a at 1300 °C will init- 
ially crystallize diopside. A line drawn from the 
CaMgSi 2 O fi apex to a, if extended beyond a, indicates 
the changes in melt composition caused by diopside 
crystallization. If the temperature continues to fall, the 
melt composition will eventually reach the boundary 
between the diopside and plagioclase ss fields (at point 
b), and here crystals of plagioclase begin to crystallize 
together with diopside. The boundary indicates the 
restricted series of melt compositions that can coexist 
with both diopside and plagioclase at the temperatures 
shown. 

To work out the composition of the plagioclase that 
crystallizes from melt b requires the use of tie-lines. But 
one must remember that tie-lines are isothermal lines 
(since two phases in equilibrium must have the same 
temperature), and for this purpose it is appropriate to 
use a second type of diagram derived from the three- 
dimensional model. This is the isothermal section shown 
in Figure 2.9d. One can visualize this section as a horiz- 
ontal slice through Figure 2.9a at a specified tempera- 
ture (see Figure 2.9c). In principle an isothermal section 
can be drawn for any temperature for which phase 



(a) (b) 




(c) (d) 




Figure 2.9 Various ways of representing crystallization in the pseudo-ternary system CaMgSi 2 0 6 -NaAlSi 3 0 8 -CaAl 2 Si 2 0 8 . (a) 
A perspective sketch of the liquidus surface in three dimensions. Elevation and contours represent temperature, (b) A plan 
view of the liquidus surface, with topography shown by temperature contours (graduated in °C). (c) 3D sketch indicating the 
construction of the isothermal section shown in (d). (d) Isothermal section at 1220 °C. Tie-lines across the two-phase fields 
link coexisting phase compositions. The arrow (a tangent to the cotectic curve in Figure 2.9b) indicates the direction in which 
the melt composition b will evolve with further crystallization. This direction is controlled by the proportion in which 
diopside and plagioclase c crystallize (point d). 
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equilibrium data are available: this one relates to 
1220 °C, the temperature of the liquidus at point b in 
Figure 2.9b. The area where the liquidus lies below the 
temperature of the section (shaded in Figure 2.9d) is a 
one-phase field where only melt is stable at 1220°C. 
The rest of the diagram can be regarded as the result of 
slicing the top off the solid model at this temperature 
(Figure 2.9c), revealing three 2-phase fields, each tra- 
versed by a family of tie-lines (summarizing the results 
of phase equilibrium experiments at 1220 °C). The 
composition of the plagioclase ss in equilibrium with 
the melt b at this temperature can be read from the dia- 
gram by following the tie-line from b to the NaAlSi 3 O s - 
Ca Al ,Si,0 K edge of the diagram (point c). 

Tie-line b-c forms one boundary of a three-phase field 
representing equilibrium between melt b, plagioclase c 
and diopside (composition CaMgSi,0 6 ) at this temper- 
ature. Any point lying within this field signifies a 
physical mixture of these coexisting phases, the pro- 
portions of which could be worked out using the Lever 
Rule. Because all possible mixtures of diopside and 
plagioclase c lie to the right of melt b (along the line 
Di-c), the crystallization of these two minerals with 
cooling causes the melt composition to migrate left- 
wards, along the boundary - the cotectic - shown in 
Figure 2.9b. This direction is indicated by the arrow in 
Figure 2.9c. Point d, collinear with the arrow, indicates 
the proportion in which diopside and plagioclase (c) 
crystallize from the melt b (Exercise 5). 

More detailed interpretation of such diagrams lies 
beyond the scope of this book. Further information can 
be found in the books by Morse (1980), Winter (2009) 
and Gill (2010). 

Review 



The great diversity of reactions and assemblages 
recorded in natural igneous and metamorphic rocks 
provides many avenues for investigating the condi- 
tions under which the rocks were formed. We have 
seen that, to analyse what such mineral assemblages 
mean in terms of pressure and temperature of forma- 
tion, we can draw on two sorts of published experi- 
mental information. The primary source is the literature 
of experimental petrology, in which one can usually track 
down a number of phase diagrams relevant to the 
assemblages in a particular rock suite. Such diagrams, 
derived from well-established laboratory procedures 



(Box 2.1) that are familiar to most petrologists, present 
phase equilibrium data in an easily understood form. 
But many diagrams of this kind refer to experiments 
on simple laboratory analogues rather than on the 
rocks themselves. The melts in Figure 2.9, for example, 
fall short of having true basaltic compositions, owing 
to the absence, among other things, of the important 
element iron. (One of the many consequences of this 
defect is that equilibria in Figure 2.9b are shifted to 
higher temperatures than would be found in a real 
iron-bearing basalt.) Thus diagrams like Figure 2.9, 
although invaluable for analysing general principles of 
phase equilibrium, do not reflect in quantitative detail 
the behaviour of more complex natural magmas and 
rocks. In some circumstances it can be helpful to carry 
out experiments on natural rock powders (see, for 
example. Gill, 2010, Figure 3.9) or comparable syn- 
thetic preparations, but the results cannot be directly 
displayed in simple phase diagrams and have less gen- 
eral application. 

There are also useful applications in petrology for 
thermodynamic data (molar enthalpies, entropies and 
volumes). Using the Clapeyron equation and Le 
Chatelier 's principle, we can predict certain features of 
phase diagrams without recourse to petrological 
experiment. Molar enthalpies and entropies of pure 
minerals are measured primarily by a completely dif- 
ferent technique called calorimetry, involving the very 
accurate measurement of the heat evolved when a 
mineral is formed from its constituent elements or 
oxides. Such methods and data are less familiar to 
most geologists, and their successful application in 
solving petrological problems requires a command of 
thermodynamic theory beyond the scope of this book. 
Thermodynamics has, however, become one of the 
most versatile tools of metamorphic petrologists, ena- 
bling them to apply experimental data from simple 
synthetic systems to complex natural assemblages. 

Further reading 



Barker, A.J. (1998) An Introduction to Metamorphic Textures and 
Microstructures. Abingdon: Routledge. 

Best, M.G. (2002) Igneous and Metamorphic Petrology. Oxford: 
Wiley-Blackwell. 

Gill, R. (2010) Igneous Rocks and Processes - a Practical Guide. 
Chichester: Wiley-Blackwell. 

Morse, S.A. (1980) Basalts and Phase Diagrams. New York: 
Springer- Verlag . 
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Exercises 



2.1 Applying the Phase Rule to the reaction shown in 
Equation 2.5, discuss the variance of points X and 
Y in Figure 2.1.1a (Box 2.1). Identify the degrees of 
freedom operating at point X. 

2.2 Why does ice floating on water tell us that the 
melting temperature of ice will be depressed at 
high pressures? 

2.3 At atmospheric pressure (10 5 Pa), the following 
reaction occurs at 520 °C: 

Ca 3 Al 2 Si 3 0 12 + Si0 2 — > CaAl 2 Si 2 0 8 +2CaSi0 3 

grossular quartz anorthite wollastonite 

( a garnet ) ( plagioclase ) 



Use the data below to plot a correctly labelled P-T 
diagram for pressures up to 10 9 Pa. 





Entropy S 
J K _1 mot 1 


Volume V 
10~ e m 3 mol~ 1 


Grossular (Ca 3 ALSi,0.,) 


241.4 


125.3 


Quartz (Si0 2 ) 


41.5 


22.7 


Anorthite (CaAI 2 Si 2 O g ) 


202.7 


100.8 


Wollastonite (CaSi0 3 ) 


82.0 


39.9 



2.4 Refer to Figure 2.5. Calculate the relative propor- 
tions of melt and crystals produced by cooling a 
melt of composition m to (a) 1400 °C, (b) 1300 °C, 
and (c) 1230 °C. What are the compositions of melt 
and plagioclase at these temperatures? (Assume 
that equilibrium is maintained throughout.) 

2.5 Plot the following rock composition in the ternary 
system CaAl 2 Si 3 0 8 -CaMgSi 2 0 6 -Mg,Si0 4 (see 



Figure 2.8): plagioclase 42.5%, diopside 25.5%, 
nepheline (NaAlSiOJ 15.0%, olivine 17.0%. (Note 
that nepheline does not appear in Figure 2.8.) What 
mineral would crystallize first from a melt corre- 
sponding to this composition? 

2.6 Refer to Figure 2.9 and its caption. In what propor- 
tions must diopside and plagioclase crystallize 
from melt b to drive its composition along the 
cotectic curve (the arrow in Figure 2.9d)? 

Calculate the compositions and proportions of 
the phases present in a solid mixture of composi- 
tion a (Figure 2.9b). What would be the equilib- 
rium assemblage for this mixture at 1220 °C? What 
would be the compositions and relative proportions 
of the phases present? 

2.7 Figure 2.10 below shows the ternary eutectic in the 
system CaMgSi,0 6 (Di)-Mg 2 Si0 4 (Fo)-Mg,Si 2 0 6 
(En) at a pressure of 20 x 10 8 Pa. To what depth in 
the mantle does this pressure correspond? What 
are the compositions (expressed as Di Fo. Ery 
where x, y and z are mass percentages) of the first 
melts to form as (a) the mixture M and (b) the 
mixture N are heated through the solidus? 



CaMgSi 2 0 6 (Di) 




(Fo) (En) 



^ 


Direction of decreasing 
temperature 


E 


Ternary eutectic 



Figure 2.10 Simplified phase relations in the system 
CaMgSi,0 6 (Di)-Mg 0 SiO 4 (Fo)-Mg 2 Si 2 0 6 (En) at a pressure 
of 20 x 10 8 Pa (for Exercise 2.7). 



Q KINETICS OF EARTH PROCESSES 



In reading Chapter 2, it is easy to fall into the trap of 
supposing that reactions between minerals always 
proceed rapidly enough to reach chemical equilib- 
rium during the time available, however short that 
may be. A little thought, however, indicates that this 
cannot be so. We have seen in Figure 1.3, for example, 
that aragonite is a metastable mineral at atmospheric 
pressure, and its survival in some outcrops of meta- 
morphic rock crystallized at high pressure is a sign of 
disequilibrium: the rate of inversion to calcite has been 
too slow for the reaction to be completed before the 
processes of uplift and erosion exposed the rock at the 
surface. The same reasoning applies to occurrences 
in surface outcrops of the mineral sillimanite, which 
is stable only at elevated temperatures (Figure 2.1). 
Examination of igneous or metamorphic rocks in thin 



section often brings to light petrographic evidence of 
disequilibrium, in such forms as mineral zoning and 
corona structures (Box 3.1). From these examples it is 
clear that the rate of progress of geochemical reactions, 
and the way they respond to different conditions, are 
factors we cannot ignore. 

The measurement and analysis of chemical reaction 
rates is called chemical kinetics, a science whose simpler 
geological applications are the subject of this chapter. 
Chemical kinetics provides the theoretical basis for 
understanding how (and why) reaction rates depend 
upon temperature, a matter of fundamental geological 
importance in view of the high temperatures at which 
magmas and metamorphic rocks crystallize. It also 
provides the algebra upon which radiometric (iso- 
topic) dating methods are based (Box 3.2). 
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Box 3.1 Disequilibrium textures 

Mineral reactions that are unable to proceed to comple- 
tion leave a rock in a state of chemical disequilibrum 
which, on the scale of a thin section, may be indicated by 
a variety of disequilibrium textures. 

Coronas 

Figure 3.1.1 shows a garnet crystal in a metamorphic rock 
(a metagabbro) which, as conditions changed during meta- 
morphism, became unstable and began to react with adja- 
cent quartz crystals to form a mixture of new minerals: 
plagioclase, magnetite (finely dispersed in the plagioclase) 
and pyroxene. The products of this metamorphic reaction 
are located at the periphery of the garnet crystal, in a 
zonal arrangement known as a corona structure. A corona 
records a mineral reaction that took place too slowly for it 
to proceed to completion (i.e. complete elimination of the 
garnet crystal) before changing conditions brought it to a 
halt. The rates of such reactions are controlled by diffu- 
sion rates and are therefore strongly temperature-depend- 
ent. What remains is a frozen-in disequilibrium texture. 

It is often impossible (as here) to represent a corona 
reaction as a balanced chemical equation between the 
minerals observed, owing to the involvement of a fluid (or 
in some cases a melt) phase which can introduce and 
remove soluble reaction components without leaving any 
visible trace in Figure 3.1.1. 




Figure 3.1.1 Corona structure rimming a six-sided garnet 
crystal in a metamorphosed gabbro. Opx and cpx stand 
for ortho- and clinopyroxene respectively. Width of field 
of view 4.5 mm. (Source: Carlson & Johnson, 1991, 
reproduced with permission of the Mineralogical Society 
of America.) 



Reaction rims 

Similar textures called reaction rims arise in various igne- 
ous rocks, due to reaction between early-formed crystals 
with later, more evolved melts. For example, olivine phe- 
nocrysts may become mantled by orthopyroxene owing to 
reaction with evolved melt richer in Si0 2 (Plate 1; see 



Defining the rate of a reaction 



It is easy to accept that some reactions proceed faster 
than others, but less easy to see how such differences 
can be expressed quantitatively. What precisely do we 
mean by the rate of a reaction? 

Consider a simple chemical reaction, for example 
that between nitric oxide (NO) and ozone (a form of 
oxygen molecule comprising three oxygen atoms, 0 3 ), 
two gaseous pollutants that occur in the troposphere 
as a result of the burning of fossil fuels. These reactants 
react with each other in equal molecular proportions to 
form the products nitrogen dioxide (NO 0 ) and ordi- 
nary oxygen (0 2 ), which are also gases: 

NO+ 0 3 — > N0 2 + 0 2 -i\ 

gas gas gas gas ' * ' 



Imagine a laboratory experiment in which gaseous 
NO and 0 3 are reacted together in a sealed vessel 
equipped with sensors that monitor the changing con- 
centrations of NO, 0 3 , N0 2 and 0 2 in the reaction ves- 
sel as the reaction progresses. (How these sensors 
work need not concern us.) The reaction consumes 
NO and 0 3 , whose concentrations (c NO and c G , each 
expressed in mol dim 3 ) therefore decrease with time as 
shown in Figure 3.1. The concentrations of the prod- 
ucts increase correspondingly as the reaction proceeds. 
The rate of the reaction at any stage is the gradient 
of the right-hand graph in Figure 3.1 at the moment 
concerned. Borrowing the symbolism of calculus 
(Appendix A): 



df 



df 



(3.2) 
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A 



also Box 2.6), or pyroxene may be rimmed by amphibole 
owing to reaction with hydrous late-stage melts at lower 
temperatures. 

Zoning 

In a rock that achieved complete chemical equilibrium 
between its phases at a given temperature, all mineral crys- 
tals would be homogeneous in composition. Igneous and 
metamorphic minerals are, however, quite commonly zoned. 
Zoning indicates that intracrystalline diffusion has failed to 
keep pace with changing external circumstances. Zoning in 
igneous minerals (see Plate 2) often reflects chemical 
evolution of the melt, with which only the rim of the growing 
crystal has maintained equilibrium (Figure 3.1.2). Zoning in 
igneous and metamorphic rocks may also be a response to 
changing physical conditions (P, T, etc.). 

Exsolution 

Perthites (Figure 2.7) and similar textures in pyrox- 
enes (Plate 3) represent the solid-state decomposition 
of a homogeneous crystal into two immiscible phases 
(Figure 2.6). Intra-crystalline exsolution lamellae like 
these have a very large area of interface with the host 
crystal. The mismatch of structure across this interface 
generates a large positive interfacial energy, a situation 
undoubtedly less stable than equilibrium segregation into 




Figure 3.1.2 Formation of a zoned crystal of plagioclase 
through cooling and compositional evolution of melt. 

separate crystals. The persistence of exsolution lamellae 
indicates that diffusion through the crystal was too slow 
to allow equilibration. 

An example of how exsolution behaviour can be used to 
measure cooling rates is given in Box 3.5. 



Figure 3.1 Composition-time curves for 
reaction 3.1. 



o 

o 



o 

o 




Because one mole of NO., is produced for every mole of 
NO that is consumed by the reaction, the gradients of the 
left- and right-hand graphs in Figure 3.1 differ only in 
their sign (negative and positive respectively). This dif- 
ference is represented by the minus sign in reaction 3.2. 



Rate equation 

If we repeated the experiment with double the concen- 
tration of NO present ( c 0s being initially the same as 
before), we would find that the initial rate is doubled. 
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Box 3.2 Kinetics of radioactive decay I: the Rb-Sr system 



The radioactive isotope of rubidium, 87 Rb (Box 10.1), 
decays to the strontium isotope 87 Sr. The kinetics of this 
nuclear reaction can be treated in the same way as a 
chemical reaction: 

87 Rb — » 87 Sr+jCT + p {321) 

‘parent’ ‘daughter’ 

isotope isotope 

where the /3-particle ([)-) and the antineutrino ( v ) are 
released by the reaction. 

The decay rate of any radioisotope is proportional to the 
number of the radioisotope (‘parent’) nuclei present in the 
sample (n ) at the moment in question. This can be written 
as a rate equation: 

do 

decay rate = - = Xn p (3.2.2) 



where n° is the number of parent nuclei initially present 
(when t = 0). The decay profile for a short-lived isotope is 
shown in Figure 3.2.1a. If we transform both sides of 
Equation 3.2.3 into natural logarithms and rearrange 
things a bit, we get: 



Ini ^ 



= Xt 



(3.2.4) 



This equation is more useful because it is linear in relation 
to time (Figure 3.2.1b). For radioisotopes that decay rap- 
idly, the decay constant can be determined by measuring 
the gradient of this graph (see Exercise 3.1). 

The half-life t, of a radioisotope is the time it takes for 

2 

n p to decay to half of its original value (n p = § n°). Thus 

ln(2/l) = Xti 

2 



Because there is only one concentration term (n ) on the 
right-hand side (unlike Equation 3.3), this is called a first- 
order reaction. X is the rate constant analogous to k in 
Equation 3.3, but in this context it is called the decay con- 
stant (same algebra, different names). 

The rate equation can be integrated to show how n p 
varies with time: 

n p = n p ° e- u (3.2.3) 



therefore t 4 = 0.6931/A, (3.2.5) 

2 

Because the decay of 87 Rb is very slow, the numerical 
value of X is extremely small: 1.42 x 1CF 11 year 1 . During 
one year only about 14 out of every million million 87 Rb 
nuclei are likely to decay. The 87 Rb remaining in the Earth 
today has survived from element-forming processes 
(Chapter 11) that took place before the solar system was 
formed 4.6 billion years ago. 



(a) (b) 




Figure 3.2.1 (a) The falling number n p of radionuclei in a sample plotted against time, (b) The function ln[n°/n ] plots as 
a straight line against time; the gradient equals X. 

V J 
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If we double the initial concentrations of both NO and 
O a , we find the initial rate is quadrupled. This suggests 
that the rate is related to reactant concentrations: 

rate = _ ~rT = k ' c No-c° 3 (3.3) 

df 

Equation 3.3 is called the rate equation for this reaction. 
Because it contains two concentration terms (c NO and 
Cq ), the reaction is said to have second-order kinetics. 
The constant k, whose numerical value is specific to 
this reaction (and to the temperature at which the 
experiment is run), is called the rate constant. The equa- 
tion predicts that as the reactants are used up the rate 
will decline, which is consistent with the flattening out 
of the slopes in Figure 3.1. 

The process of radioactive decay can be analysed in 
a similar manner (Boxes 3.2 and 3.3). 

Heterogeneous reactions 

Reactions like 3.1 that take place within a single 
phase (in this case a homogeneous gas mixture) are 
called homogeneous reactions. Nearly all reactions 
of geological significance, on the other hand, are het- 
erogeneous reactions, involving the participation of 
two or more phases (minerals, melts, solutions ...). 
Because they require the migration of components 
across the interface dividing one phase from another, 
the formulation of rate equations for heterogeneous 
reactions is much more complicated than for homo- 
geneous reactions. 

The most obvious consequence of involving two 
phases in a reaction is that the surface area of their 
interface becomes a variable in the rate equation. 
Interfacial surface area is determined chiefly by parti- 
cle size. The surface area of a cube 1 cm across is 6 cm 2 
(six sides each of 1 cm 2 area). Cutting the cube in half in 
each direction produces eight cubes, each 0.5 cm across 
and each having a surface area of 6 x 0.5 2 = 1 .5 cm 2 . The 
total volume of all the cubes together is unchanged 
(lcm 3 ) but the total surface area has increased from 
6cm 2 to 8 x 1.5 = 12cm 2 . Dividing the original cube into 
1000 cubelets each of 0.1cm size would increase the 
total area to 60 cm 2 , while reducing to particle sizes 
equivalent to silt and clay sediments would increase 
their surface area to 3000 and 60,000 cm 2 respectively. 
Particle or crystal size, because it determines the area 
of contact between phases, has a profound effect on the 



rate of a heterogeneous reaction. This is one reason 
why diesel fuel injected as a fine spray into an engine 
reacts explosively with air, whereas the bulk liquid 
burns much more slowly. 

The condition of the interface is also a very import- 
ant factor. The rate of inversion of aragonite to calcite, 
for example, is greatly accelerated by the presence of 
traces of water along the grain boundaries. Surface 
chemistry has many important applications in the 
chemical industry and in mineral processing (for 
example, the use of a frothing agent to optimize the 
separation of ore minerals by flotation). 

Mechanical factors come into play as well. When 
a solid dissolves in still water, the aqueous phase 
surrounding it becomes locally saturated, and this 
impedes further solution until diffusion has distrib- 
uted the dissolved species more evenly. Dissolution of 
sugar in coffee can therefore be accelerated by the use 
of a teaspoon to promote homogenization, and natural 
forms of agitation can be correspondingly effective in 
the marine environment. Experiments show that the 
rate at which calcite dissolves in water can be rep- 
resented like this: 

Rate = Mod |fC° - (c Ca2+ (c ccv _ )"| (3.4) 

The c terms refer to concentrations of ions in solu- 
tion, K° and k are constants, A is the total surface area 
of the calcite phase present, and a is the experimental 
stirring rate (which appears as the cube root for reas- 
ons that need not concern us). No doubt the effect of 
natural wave-agitation is still more complicated. This 
equation illustrates how rapidly the complexities mul- 
tiply when even the simplest heterogeneous reactions 
are studied kinetically. 

Temperature-dependence of reaction rate 



Everyday experience tells us that chemical reactions, 
whether homogeneous or heterogeneous, speed up as 
the temperature is raised. Epoxy adhesives cure more 
quickly in a warm oven. Conversely, the very fact that 
we use refrigerators and freezers to preserve food 
indicates that biochemical reactions slow down at 
lower temperatures. Quantitatively the temperature 
effect which these examples illustrate is quite pro- 
nounced: many laboratory reactions roughly double 
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Box 3.3 Kinetics of radioactive decay II: the U-Th-Pb system 



Each of the naturally occurring isotopes of the trace ele- 
ments uranium ( 235 U and 238 U) and thorium ( 232 Th) decays 
through a complex series of intermediate radioactive 
nuclides to an isotope of lead (Pb). This is illustrated for 
the decay of 238 U to 206 Pb in Figure 3.3.1; similar flow dia- 
grams can be drawn for the decay of 235 U to 207 Pb and for 
232 Th to 208 Pb. Together these decay chains form the basis 
of U-Th-Pb radiometric dating. 

Despite its complexity, the overall decay process in 
Figure 3.3.1 conforms to first-order kinetics, because the 
first step in the process (to 234 Th, a short-lived radioactive 
isotope of thorium) happens to be the slowest. The kinetic 
complexities of the subsequent branching decay series are 
immaterial because the rate of the whole process is con- 
trolled by this one rate-determining step, just as the flow of 



water from the end of a hose can be controlled by adjusting 
the tap supplying it. This phenomenon is not limited to radi- 
oactive decay: the kinetics of some complex chemical reac- 
tions are also controlled by a slow, rate-determining step. 

For every uranium or thorium nucleus that decays to 
lead within the Earth, between 6 and 8 alpha particles are 
released. By capturing electrons the a-particles become 
4 He atoms, which form the bulk of the helium flux escaping 
from the Earth’s interior. 

All but one of the nuclides involved in Figure 3.3.1 are 
radioactive solids, likely to be retained within the mineral 
hosting the original U. The sole exception is 222 Rn, one of 
the isotopes of the inert gas radon, whose mobility pre- 
sents an environmental hazard in areas underlain by high- 
U-Th rocks like granites, as discussed in Box 9.9. 




Figure 3.3.1 The chain of radioactive decay steps by which 238 U (also written ‘uranium-238’) decays to 206 Pb (lead-206). 

V J 



their reaction rates when the temperature is raised by 
just 10°C (see Exercise 3.2 at the end of this chapter). 
The temperature-dependence of reaction rates is par- 
ticularly significant for geological processes, whose 
environments can vary in temperature over many 
hundreds of degrees. 

Most reaction rates vary with temperature in the 
manner shown in Figure 3.2a. The Swedish physical 



chemist Svante August Arrhenius 1 showed in the 
late 1880s that this behaviour could be represented 



1 Arrhenius is also notable for being the first scientist to postu- 
late (in 1896) a climatic 'greenhouse effect' arising from the 
presence of CO, in the atmosphere, and for recognizing even 
in the nineteenth century that mankind's burning of fossil car- 
bon contributed to global warming. 
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Figure 3.2 Variation of the rate constant with temperature, (a) k plotted directly against T (in kelvins). (b) k plotted against 
reciprocal temperature. T _1 has been multiplied by 10 3 to give a more convenient number scale. The reciprocal temperature 
scale is graduated directly in kelvins along the top margin. (This axis shows the recognized way of writing 10 3 /T with T 
expressed in kelvins). 



algebraically by expressing the rate constant k for a 
reaction in terms of an exponential equation: 

k = Ae~ E ‘ /RT (3.5) 

This has become known as the Arrhenius equation and it 
has a number of important applications in geology. R is 
the gas constant (8.314 J mol -1 K -1 ) and T represents 
temperature in kelvins. A and E a are constants for the 
reaction to which the rate constant refers (varying from 
one reaction to another). A is called the pre-exponen- 
tial factor, and it has the same units as the rate constant 
(which depend on the order of the reaction concerned). 
The constant E a is called the activation energy of the 
reaction and has the units J mol -1 . 

Chapter 1 showed that physical and chemical pro- 
cesses often encounter an energy 'hurdle' that impedes 
progress from the initial, high-energy (less stable) state 
to the lower-energy configuration in which the system is 
stable (Figure 1.4). The hurdle is illustrated for a hypo- 
thetical chemical reaction in Figure 3.3. It arises because 
the only pathway leading from reactants to products 
necessarily passes through a higher-energy transition 
state. For a reaction involving the rupture of one bond 
and the establishment of a new one, this involves the for- 
mation of a less stable intermediate molecular species 




Progress of reaction *. 



Figure 3.3 Energetics of a hypothetical reaction AB + C -> 
A+BC. The vertical axis represents the free energy of the 
system. E a is the activation energy, which is released again 
(except in the case of endothermic reaction, when only a 
part is released) on completion of the molecular reaction. 

(the activated complex in Figure 3.3). The activation energy 
E a in the Arrhenius equation may be visualized as 
the 'height' of the hurdle (in free energy units) relative to 
the initial reactant assemblage. Box 3.4 explains how the 
activation energy can be understood on the atomic scale. 
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Box 3.4 What activation energy means on the atomic scale 



In a chemical reaction (Figure 3.3): 

AB + C -> A + BC 

the established A-B bond must be weakened (stretched) 
before a new bond (B-C) can begin to form. The energetics 
of the A-B bond are shown in Figure 3.4.1a. The process 
begins with AB in its most stable configuration (inter- 
nuclear distance = r A °_ B ). Energy is required to stretch the 
A-B bond to the stage when formation of the B-C bond 
becomes an equally probably outcome (i.e. sufficient to 
form the activated complex A - B - C); this energy input 



constitutes the activation energy E a (Figures 3.3 and 
3.4.1c). The whole reaction from AB to BC can be visual- 
ized by considering the energy-distance curves of both 
molecules, Figures 3.4.1a and b, ‘back to back’ as in 
Figure 3.4.1c. Note that the AB bond need not be com- 
pletely broken before the BC molecule can begin to form. 

The explanation of the activation energy in reactions 
between ionic compounds is slightly different, but is still 
associated with the need to disrupt one arrangement of 
atoms or ions before another more stable arrangement 
can be adopted. 




Figure 3.4.1 (a) and (b). Bond energy versus nuclear separation r A B and r A B for molecules AB and BC. (c) Bond energy 
and activation complex in the transition from AB + C to A + BC. 
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Figure 3.4 Molecular kinetic-energy distribution for two 
temperatures. The shaded areas show the portions of each 
distribution that lie above a specified energy E x . This 
proportion is greater at the higher temperature. 

For a collision between reactant molecules to lead to 
the formation of product molecules, it must derive suf- 
ficient thermal energy from the participants to gener- 
ate the activated complex. From theoretical calculations 
one can predict the kinetic energy distribution among 
reactant molecules at a given temperature T: the results 
are illustrated for two alternative temperatures in 
Figure 3.4. One can show that the proportion of molec- 
ular encounters involving kinetic energies greater than 
some critical threshold energy E x (shaded areas in 
Figure 3.4) is given by: 

Proportion of molecular collisions /q s\ 

exceeding E x x e~ E * /RT 

The term on the right is called the Boltzmann factor. It 
appears in the Arrhenius equation as a measure of the 
proportion of the reactant molecule collisions that pos- 
sess sufficient energy (i.e. greater than E ) at the tem- 
perature T to reach the transition state and thereby to 
complete the reaction. The form of the Boltzmann fac- 
tor shows that an increase of temperature will shift 
the energy distribution towards higher energies 
(Figure 3.4), so that a greater proportion of reactant 
molecules can collide with energies exceeding E and 
surmount the energy hurdle, like water flowing over a 
weir. In other words the reaction rate will increase. But 
reducing the temperature will inhibit the reaction, 
because -E /RT becomes a larger negative number 



and makes the Boltzmann factor and therefore k 
smaller. 

The value of the activation energy tells us how sensi- 
tive the rate of a reaction will be to changes in temper- 
ature, and measuring this temperature-dependence 
provides the means by which E a can be determined 
experimentally. (Because it is not a net energy change 
between reactants and products like AH, E a cannot be 
measured calorimetrically.) Transforming both sides of 
Equation 3.5 into natural logarithms ('In' - see 
Appendix A), one gets: 

lnfc = - — ■ — + lnA (3.7) 

R T 

which has a linear form: y = m . x + c (see Appendix A). 
Thus if In k {'y') is plotted against 1/T ('x , T being in 
kelvins), a straight line should be observed as shown 
in Figure 3.2b and discussed in Appendix A. The 
gradient of this line is -EJR. The numerical value 
of the activation energy can therefore be established 
by repeating the rate experiment at a number of pre- 
determined temperatures, and plotting the rate con- 
stants obtained in the form shown in Figure 3.2b. Such 
a diagram is called an Arrhenius plot. An Arrhenius plot 
for reaction 3.1 is illustrated in Figure 3.5a. A geochem- 
ical example is shown in Figure 3.5b. See also Exercise 
3.2 at the end of the chapter. 

Photochemical reactions 

Is a rise in temperature the only mechanism by which 
reactant molecules can surmount the activation energy 
barrier? A number of gas reactions known to take place 
in the stratosphere show that energetic photons from 
the Sun provide an alternative means of energizing 
molecules into reacting. One example is the formation 
of stratospheric ozone (0 3 ): 

O, UVphoton > 0+0 (3.8) 

The intense flux of solar ultraviolet (UV) radiation - 
see Box 6.3 - in the upper atmosphere causes a small 
proportion of oxygen molecules to disintegrate into 
separate oxygen atoms which, by virtue of their dis- 
rupted chemical bonds (represented by the symbol 
0*),are highly reactive: they effectively sit on top of the 
energy barrier shown in Figure 3.3. Collision of either 
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(b) 




Figure 3.5 (a) Arrhenius plot for reaction 3.1. The rate constant k is expressed in dm 3 mol _1 s _1 . (b) Arrhenius plot for the 
photochemical oxidation of pyrite after here the rate k is expressed in units of pmol nr 2 min. Note that temperature is 
expressed in kelvins (K) in each case. (Sources: After Zhang (2008). © Princeton University Press, based on data of Borders & 
Birks (1982); Schoonen et al. (2000). Reproduced with the authors' kind permission.) 



of these oxygen free radicals with another oxygen 
molecule generates an ozone (tri-oxygen) molecule, as 
shown in reaction 3.9. 

0-+0 2 ^0 3 (3.9) 

Chemical reactions like 3.8 that are initiated by 
energetic photons are described as photodissocia- 
tion reactions. In kinetic terms they differ from con- 
ventional reactions in that their rate constants depend 
not on reactant temperature (as in Equation 3.5) but 
on the UV photon flux to which the reactants are 
exposed. 

Solar UV radiation in a different wavelength range 
stimulates a reverse reaction: 

a, second UV photon > + Q . (3.10) 

This is the key stratospheric reaction that strongly 
attenuates the incoming UV flux in this wavelength 
range and protects living things on the Earth's surface 
from the most damaging effects of solar UV radiation. 
This filter is highly efficient: we often refer to the 
'ozone layer' protecting us, yet remarkably (in view of 
the vital protection all living things on Earth derive 
from it) the steady-state ozone concentration in the 
stratosphere rarely exceeds 5 ppm. 



Ozone also forms photochemically in the lower trop- 
osphere, but here the intensity of solar UV is much 
lower, and ozone formation can only occur when catal- 
ysed by nitrogen oxides arising from traffic pollution 
(a factor in the infamous 'photochemical smog' that 
afflicts many of the world's sunniest cities today, espe- 
cially those at higher altitude like Mexico City). 



Diffusion 



Any process that requires the input of thermal energy 
to surmount an energy hurdle - a thermally activated 
process - will show an Arrhenius-type temperature 
dependence (Equation 3.5 and Figure 3.2). This prop- 
erty is characteristic of a number of physical processes 
in geology in addition to geochemical reactions. 

When a component is unevenly distributed in a 
phase - whether solid, liquid or gas - so that its con- 
centration in one part is higher than in another, ran- 
dom atomic motions tend to 'even out' the irregularities 
over time, leading to a net migration, called diffusion, 
of the component 'down' the gradient to regions of 
lower concentration. Given sufficient time, diffusion 
leads eventually to a homogeneous distribution 
(labeled t^ in Figure 3.6a). In fact, atoms diffuse 
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Figure 3.6 (a) The diffusion of ions in response to a non-uniform concentration distribution, (b) Arrhenius plot of diffusion 
coefficients for the diffusion of cobalt ions in silicate melts. The activation energies measured from the slopes of the two graphs are: 

basalt: 220 kj mob 1 

andesite: 280 kj mob 1 

The units of D are m 2 sA Temperature is expressed in kelvins (K). (Source: Adapted from Lowry et at. (1982). 

Reproduced with permission from Springer Science and Business Media.) 



randomly even through a homogeneous substance, 
but only where a concentration gradient exists is a net 
flow of chemical components observed. 

Imagine a crystal with an abrupt internal discont- 
inuity in the concentration of a component i, as shown 
in Figure 3.6a. Suppose the crystal is maintained at a 
constant temperature T, high enough for solid-state 
diffusion to occur at a significant rate. If we were to 
measure the distribution of component i on several 
successive occasions t v t 2 , t y we would see the devel- 
opment of a progressively smoother concentration 
profile, leading eventually to a uniform distribution of 
i (at f in Figure 3.6a). These changes point to a net flux 
of component i from left to right through the plane of 
the original discontinuity. The magnitude of the flux in 
moles per second will depend on the surface area of 
this interface. One therefore expresses the flu xf as the 
amount of component i (in moles) that migrates 
through a unit area of the plane per second, so the 
units are 'moles per square metre per second' 
(mol m 2 s '). Common sense suggests that the flux will 
depend upon the steepness of the concentration grad- 

An 

ient l ; so no net flux will occur if the concentration is 

dx 



uniform and the gradient zero, whereas a high flux 
will occur if the gradient is steep. The 19th century 
German physicist Adolf Fick backed up this hunch 
theoretically in 1855, showing that: 

fi=~D,p- (3.11) 

dx 

This equation is known as Fick's First Law of Diffusion. 
The negative sign indicates that the direction of the net 
flu xf is down the concentration gradient (i.e. towards 
the right in Figure 3.6a). The constant D ; is called the 
diffusion coefficient for the species i in the crystal con- 
cerned (at a given temperature). The units of concen- 
tration c are 'moles per cubic metre' (mol nr 3 ), so the 
concentration gradient dc / dx will have units of moles 
per cubic metre per metre in the x direction ([mol nr 3 ] 
nr 1 = mol nr 4 ) from which it is simple to show that the 
units of D must be rrrs 1 (see Exercise 3.5). 

Many experimental determinations of diffusion 
coefficients have been made for various elements in a 
range of silicate materials at various temperatures. 
Figure 3.6b shows how the measured diffusion 
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coefficient for cobalt (D ) changes with temperature 
and composition in basalt and andesite melts. The data 
are presented in the form of an Arrhenius plot, in 
which they lie on straight lines. Thus one can express 
the temperature dependence of D in terms of an equa- 
tion similar to Equation 3.5: 

D = D 0 e~ E ‘ /RT (3.12) 

or in log form: 

InD = — h In D„ (3.13) 

R T o v > 

These equations are identical to Equations 3.5 and 3.7, 
except that the pre-exponential factor is written D g . 
Thus, although we think of diffusion as a physical 
phenomenon, it resembles a chemical reaction in being 
governed by an activation energy E a . Like a person 
caught up in a dense crowd, the diffusing atom or ion 
must jostle and squeeze its way through the voids of 
the melt (Chapter 9), and pushing through from one 
structural site to the next presents an energy hurdle 
which only the more energetic ('hotter') atoms can 
surmount. 

Measured diffusion coefficients differ from one 
element to another (Co as opposed to Cr, for example) 
for diffusion in the same material at the same temper- 
ature. Figure 3.6b shows that the diffusion coefficient 
for one element in a melt will also vary with the com- 
position of the melt (Box 9.2). 

Solid-state diffusion 

From the point of view of diffusion, a crystalline solid 
differs from a melt in several important respects. 
Atoms diffusing through a silicate melt encounter a 
continuous, isotropic medium (D is independent of 
direction) that is a relatively disordered structure. Most 
crystalline solids, on the other hand, are polycrystal- 
line aggregates that offer two routes for diffusion: 
within and between crystals. 

Volume (intra-crystalline) diffusion 

Volume diffusion through the three-dimensional vol- 
ume of the constituent crystals is similar in general 
terms to diffusion through a melt. The main difference 
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Figure 3.7 Comparison of volume diffusion rates in silicate 
melts and olivine crystals. (Source: Data from Henderson 
1982). 

is that the crystals have more closely packed, ordered 
atomic structures than melts (Box 8.3), and diffusion 
through them is much slower (lower D). Figure 3.7 
contrasts the diffusion behaviour of similar metals in 
melts and in olivine crystals (in which the diffusion 
coefficients are lower by a factor of about 10 s ). Notice 
that diffusion in olivine takes place more readily along 
the crystallographic z-axis than along the y-axis. 
Crystallographic direction is an important factor in 
diffusion through anisotropic crystals, reflecting the 
internal architecture of the crystals (Chapter 8). 

Grain-boundary (inter-crystalline) diffusion 

Grain-boundary diffusion exploits the structural dis- 
continuity between neighbouring crystals as a channel 
for diffusion. This is much more difficult to quantify 
because the rate of diffusion depends upon: 

(a) the grain size of the rock: in a fine-grained rock, 
the total area of grain boundaries is larger in rela- 
tion to the total volume of the rock, and diffusion 
will be easier; 

(b) the microscopic characteristics of the grain bound- 
aries, e.g. the presence or absence of water. 
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Grain-boundary diffusion appears to be much more 
rapid in most circumstances than volume diffusion. 
Grain boundaries provide the main conduit through 
which volatile species penetrate the volume of a rock 
during metamorphism and hydrothermal alteration. 
Open fractures and faults do not occur in the deep crust 
owing to the high temperature and confining pressure. 
Under these conditions fluid migration is concentrated 
along shear zones, where deformation has led to a local 
reduction of grain size and often to a strong foliation, 
both of which promote diffusion and fluid movement 
because of the increase in grain-boundary area. 

Box 3.5 illustrates how diffusion studies can be 
applied to estimate the cooling rates experienced by 
iron meteorites. 



Table 3.1 Illustrative viscosity values (in Pas) for everyday 
liquids and for dry silicate melts 



Substance 


Viscosity at room 


Melt viscosity 




temperature 


above liquidus 




(-25° C) 


temperature 


Water 


0.001 




Motor oil SAE 50 


0.5 




Egg white 


2.5 




Treacle 


20 




Basalt melt (dry)* 




7.5-150 


Smooth peanut butter 


250 




Glazier's putty 


~10 5 




Rhyolite melt (dry)* 




~10 10 



*Melt viscosity increases if suspended crystals or bubbles are present. 



Melt viscosity 



The flow of silicate melts is another kind of geological 
process where the rate varies strongly with temper- 
ature. Like treacle, silicate melts become less viscous as 
the temperature is raised. 

A liquid flows in response to shear stress applied to 
it, which generates a velocity gradient in the liquid. 
The velocity of water flowing in a pipe of radius r, for 
example, increases from zero at the wall to a maximum 
v at the centre, a mean gradient of v/r. For most liquids 
the velocity gradient dv/dz is proportional to the 
applied shear stress o: 



dz; _ ct 
dz q 



(3.14) 



The parameter y in Equation 3.14 is called the viscosity 
of the liquid. Because viscosity measures resistance to 
flow, it is like an inverse rate constant: a low viscosity 
indicates a runny liquid capable of rapid flow, whereas 
a high value signifies a 'stiff' or viscous one that flows 
only slowly. 

As to units of measurement, shear stress a is meas- 
ured in Pa = N m -2 (Appendix A, Table A2), and dv/dz 
is measured in (ms _1 )m _1 = s _1 . Rearranging Equation 
3.14 gives r] = a /dv/dz, so the units in which viscosity 
is measured are Pa/(s _1 ) = Pas. To provide a feel for 
how widely viscosities vary in nature. Table 3.1 
provides some illustrative values. 



The viscosity of an erupted lava affects the style of a 
volcanic eruption, and its value depends upon both 
melt composition and temperature. Siliceous lavas like 
rhyolite have viscosities several orders of magnitude 
higher than basalts at their relevant liquidus tempera- 
tures (Table 3.1), as explained in Box 8.3. Figure 3.8 
shows both this compositional effect and the tempera- 
ture variation. By analogy with Figures 3.5 and 3.6b, 
viscosity is plotted in a logarithmic form that varies 
linearly with inverse temperature. Each of these 
straight lines can be described by an equation involv- 
ing an activation energy E a : 



1 

h 



1 C ~E JRT 

ho 



(3.15) 



or using natural logarithms : In 




- — — + In 
R T 




(3.16) 



The resemblance of Equation 3.16 to Equation 3.13 
tells us that the flow of a silicate melt - like diffu- 
sion - requires atoms or molecules to jostle past each 
other, surmounting energy hurdles as they do so. 
The collective effect of these energy hurdles is 
expressed in the flow activation energy E a , the value 
of which can be determined by measuring the gradi- 
ent of each line (= -E /R). The activation energies 
for the flow of silicate melts (which Figure 3.8 shows 
to be, like viscosities themselves, higher for siliceous 
melts than for low-silica melts) are generally in 
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Box 3.5 Diffusion and cooling rate: applications to meteorites 



At temperatures below 900°C, the continuous high-tem- 
perature solid solution between metallic iron and nickel 
divides into two solid solutions, separated by a two-phase 
region below a solvus (Figure 3.5.1a). Accordingly many 
iron meteorites exhibit complex exsolution intergrowths of 
two separate Fe-Ni alloys, kamacite and taenite, devel- 
oped during cooling through this solvus (Figure 11.1.1). 

Consider the cooling of a crystal of alloy X which at 
nearly 900°C is a homogeneous metallic solid solution. At 
about 700°C the alloy encounters the solvus, below which 
taenite becomes supersaturated with Fe, and expels it in 
the form of platelets of a separate kamacite phase, whose 
composition is indicated by drawing a tie-line to the other 
limb of the solvus (Ni ~ 6%). The upper inset shows a pro- 
file of the Ni content in a representative cross-section of 
the crystal (on the scale of a few mm) at 680°C: slender 
lamellae of low-Ni kamacite have appeared in the initially 
homogeneous taenite. 

The Fe-Ni solvus is unusual in that both limbs slope in 
the same direction (cf. Figure 2.6). As the temperature 
falls, therefore, both immiscible phases grow more Ni-rich. 



The lever rule (Box 2.3) indicates that their relative propor- 
tions must also change: the middle and lower insets show 
that kamacite lamellae grow at the expense of taenite, 
forming ever-thicker plates. 

This process relies on the diffusion rates of Fe and Ni 
atoms in the two phases. Laboratory experiments show 
that diffusion rates are slower in taenite than in kamacite, 
and thus Ni is expelled from kamacite more rapidly than it 
can diffuse into the interior of the adjacent taenite crystal. 
Consequently the Ni distribution in taenite lamellae devel- 
ops an M-shaped profile, with Ni concentrated at the 
edges. The more rapidly the meteorite is cooled, the more 
pronounced is the central dip in the Ni profile. Calculations 
based on diffusion profiles allow estimation - from the 
shape of the M-profiles in iron meteorites - of the cooling 
rates they experienced during the early development of 
the solar system (Figure 3.5.1b). These estimates, com- 
monly between 1°C and 10°C per million years, suggest 
that iron meteorites are derived from relatively small par- 
ent bodies (diameter < 400 km); large planetary bodies 
would cool more slowly (Flutchison, 1983). 




Figure 3.5.1 (a) Phase relations for Fe-Ni alloys showing the subsolidus fields of kamacite and taenite and the 2-phase 
field (ruled) between them (see Goldstein and Short, 1967). (b) Calculated diffusion profiles across a taenite platelet as 
a function of cooling rate, based broadly on Wood (1964). (Source: Adapted from Hutchison (1983)). See http://www. 
psrd.hawaii.edu/April07/irons.html 

V J 
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the same range as for diffusion (see captions to 
Figures 3.6 and 3.8). 

The deformation of crystalline rocks, although a 
much more complicated process, shows a similar 
dependence on temperature: 

strain rate = (Ae~ E,/RT (3.17) 

where a is again shear stress, and N is a constant 
depending on details of the deformation process (its 
value is typically about 8). The factor in brackets is 
analogous to inverse viscosity. It is not given that name 
because a, instead of appearing in linear form, is raised 
to the power N. Activation energies are nonetheless in 
the same range as those for viscous flow of melts. 

Persistence of metastable minerals: 
closure temperature 



An equilibrium Earth would be extremely boring (e.g. there would 
be no life, no 0 2 in the air, no plate tectonics, etc.). Disequilibrium is 
what makes the Earth so diverse and interesting. Without kinetic 
barriers there would be no geochemists to study kinetics or science, 
because all human beings, and in fact all life forms, 'should burst 
into flames'! (Zhang, 2008) 

Paradoxically, chemical equilibrium in geological pro- 
cesses would be of less interest to the petrologist were it 
not for the intervention of disequilibrium. Consider an 
argillaceous rock undergoing metamorphic recrystalli- 
zation in conditions corresponding to the kyanite-sill- 
imanite phase boundary (Figure 2.1). Once equilibrium 
is achieved, kyanite and sillimanite will coexist stably 
in the rock. If, following uplift and erosion, the rock is 
found at the surface still containing kyanite and sill- 
imanite, their coexistence and textural relationship will 
testify to the high-temperature conditions under which 
the rock crystallized and will give the petrologist an 
indication of the conditions of metamorphism. Yet the 
current state of the rock is plainly one of disequilib- 
rium, as sillimanite is now well outside its stability 
field; it remains only as a metastable relic of a former 
state of equilibrium. 

Metamorphic equilibrium seems generally to keep 
pace with changing conditions more effectively in pro- 
grade reactions (those involving a temperature increase) 
than during the waning stages of regional metamor- 
phism. One reason is that volatile constituents such as 
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Figure 3.8 Arrhenius plot of viscosities (in Ns nr 2 ) for some 
silicate melts; Viscosity is the inverse of the flow rate-constant 
so it has been plotted here in reciprocal form. Temperature 
is expressed in kelvins (K). The slope of the graphs increases 
with silica content, reflecting an increase in activation 
energy as follows: 

Basalt melt: 230 kj mol -1 

Rhyolite melt: 350 kj mol -1 

Pure silica melt: 500 kj mob 1 

(Source: data from Scarfe (1978) and other authors cited 

therein). 

water vapour, whose presence accelerates recrystalliza- 
tion, will be more abundant during the prograde stage 
owing to the dehydration reactions taking place. Once 
volatiles have been lost from the system, they are no 
longer available to facilitate the retrograde (falling-tem- 
perature) stages of metamorphism that follow. In the 
absence of volatiles it is likely that temperature asserts 
the dominating influence on metamorphic reaction rates. 

As metamorphic or igneous rocks cool, they pass 
into a temperature range of increasing kinetic paraly- 
sis, in which reaction rates fall behind changing cir- 
cumstances (as illustrated in Box 3.5). Eventually a 
temperature is reached (depending on the reaction 
being considered) at which diffusion rates become 
substantially slower that the cooling rate of the rock. 
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and reaction has effectively ground to a halt. This is 
called the blocking or closure temperature. 

The closure temperature is an important concept in 
geochronology. In the potassium-argon dating method 
(Chapter 10) the age of a rock is established by measur- 
ing the minute amount of 40 Ar (argon is a gas) that has 
accumulated in a potassium-bearing mineral from the 
decay of 40 K. The method therefore depends critically 
on the ability of the relevant mineral grains to imprison 
this intra-crystalline argon component that, being a 
gas, tends to escape. At temperatures above the closure 
temperature for argon diffusion, the 40 Ar atoms diffuse 
to grain boundaries and escape. A K-Ar age determ- 
ination on an igneous rock therefore records not the 
age of intrusion or eruption, but the date when the 
rock had cooled to temperatures low enough for the 
rate of diffusion of 40 Ar out of the grains of the potas- 
sium mineral to be insignificant. A later phase of meta- 
morphism in which the rock is again heated above the 
closure temperature would discharge the 40 Ar accumu- 
lated up to that point, and the isotopic clock would 
thereafter record a 'metamorphic age', indicating the 
time when the rock body had again cooled below the 
closure temperature. 

Review 



A chemical reaction can be visualized in free energy 
terms as a journey leading from one valley, the domain 
of the reactants, to another where the product species 
form the dominant population (Figure 3.3). The route 
from one valley to the other leads across a high pass in 
free energy space, the 'transition state', and reactant 
molecules that cannot summon sufficient energy to 
traverse this high ground will not reach their 'destina- 
tion'. The key factor is the availability of energy: if 
energy is plentiful, the traffic over the pass will be 
heavy (the reaction rate will be high); if energy is hard 
to come by, many reactant molecules will be forced to 
turn back before attaining the pass. 

What sources of energy drive geochemical reac- 
tions? For reactions taking place in the Earth's interior, 
reactants must rely on the kinetic energy possessed by 
their constituent atoms and molecules. Accordingly 
the temperature plays a crucial role in determining 
the rates of geochemical reactions, as enshrined in 
the Arrhenius equation. Reactions involving silicate 



minerals are characterized by large activation energies, 
and therefore only at high temperatures is there a sig- 
nificant population of reactant 'molecules' possessing 
the kinetic energy required to surmount the activation 
energy hurdle (Figure 3.4). Disequilibrium textures in 
many rocks (Box 3.1) testify to the rapid slo wing-down 
of chemical reactions as temperature falls, and at sur- 
face temperatures disequilibrium is the rule rather than 
the exception: the persistence of Fe 2+ minerals on the 
Earth's surface, where atmospheric oxygen makes Fe 3+ 
the stable form of iron, is one obvious example. Both 
the flow of silicate melts (Figure 3.8) and diffusion are 
also strongly temperature-dependent (Figures 3.6 and 
3.7), suggesting they too involve an activation step. 
The marked slowing-down of diffusion with falling 
temperature, and its effective cessation at the closure 
temperature, are essential requirements for radiomet- 
ric dating. 

Many reactions in the atmosphere, however, rely on 
another energy source, the Sun. Solar photons, particu- 
larly those of UV wavelengths, are energetic enough to 
tear apart chemical bonds in molecules like 0 2 , 0 3 , NO z , 
FI,0 2 , CFC1 3 and HCHO (formaldehyde). The result of 
such photodissociation reactions is often the formation 
of free radicals such as O, HO*, H*, HOp, NOp Cl*, 
CIO* and HCO*. Free radicals, possessing unpaired 
electrons (represented by the symbol '•'), are highly 
reactive: through photodissociation they have already 
reached the energy 'pass' and thus they are capable of 
initiating gas reactions in a manner unrestrained by 
considerations of activation energy. For example: 

CH 4 + HO* -» H 3 C • +H 2 0 (3.18) 

This reaction incidentally illustrates the key role of the 
hydroxyl free radical HO- as an atmospheric cleansing 
agent, removing many forms of pollution (including 
the potent greenhouse gas methane, CH 4 ) from the air 
we breathe. 



Further reading 
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Exercises 



3.1 The table below gives the Geiger counter count-rate 
at various times during an experiment on radioac- 
tive iodine-131 ( 131 I). Treating the decay of 131 I as a 
chemical reaction in which the count rate may be 
taken as being a measure of the concentration of 131 I, 
confirm that the reaction is first-order and determ- 
ine the decay constant. From the linear graph you 
have drawn, calculate the half-life of 131 I (Box 3.2). 



Time (hours) 


Counts (sec f) 


0 


18,032 


25 


16,410 


50 


15,061 


100 


12,590 


200 


8789 


300 


6144 


400 


4281 


500 


3002 



3.2 There is a rule of thumb that says that the rate 
of many room-temperature chemical reactions 
approximately doubles when the temperature is 



increased by 10°C. Calculate the activation energy 
for a reaction that exactly conforms to this relation- 
ship. (Gas constant R = 8. 3 143 JK -1 mol -1 .) 

3.3 Use the viscosity measurements given below to 
verify that the Arrhenius equation is applicable to 
the flow of a silicate melt, and determine the acti- 
vation energy. (Hint: viscosity is the resistance to 
flow, to which flow rate is inversely proportional.) 



Temperature 


Viscosity of rhyolite 


(°C) 


melt (Nsm~ 2 ) 


1325 


2042 


1345 


1585 


1374 


1097 


1405 


741 



(Gas constant as given in Exercise 3.2.) 

3.4 Calculate the half-life of 87 Rb(k 87Rb =1.42x10 "year 1 ). 
What percentage of the 87 Rb incorporated in the 
Earth 4.6 x 10 9 years ago has now decayed to 87 Sr? 

3.5 Demonstrate that the units of diffusion coefficient 
D i in Equation 3.11 are m 2 s -1 , when concentration c 
is expressed in mol m -3 and the net flux/' is expressed 
in molm -2 s -1 . 



4 AQUEOUS SOLUTIONS AND 
THE HYDROSPHERE 



The great importance of water and aqueous solutions 
on the surface of the Earth barely needs pointing out. 
Water is the principal agent of erosion and of the 
transportation of eroded materials, either by mechani- 
cal or chemical means. The world's oceans are the pri- 
mary medium for sedimentation, they act as a global 
chemical repository for many substances of geological 
significance, and they play a crucial part in moderat- 
ing the climate and supporting life on the planet (not 
to mention the role water plays in living things). Water 
also has important functions in the Earth's interior: 
ore transport, rock alteration and metamorphism all 
involve the migration of hot aqueous fluids through 
the crust. The chemistry of aqueous solutions - the 
subject of this chapter - is a vital factor in all these 
geological processes. 



A solution consists of two kinds of constituent which 
are given confusingly similar names: 

• solute refers to the dissolved species (e.g. the salt in 
a saline solution); a solution may contain several sol- 
utes, such as when sodium chloride and potassium 
nitrate are dissolved in the same solution. Seawater 
is a complex example of such a mixed solution. 

• solvent refers to the medium in which the solutes 
are dissolved. Aqueous solutions are those in which 
water is the solvent (Latin aqua : water). 

The basic jargon of solution chemistry is reviewed in 
Appendix B and key terms are defined in the Glossary. 
Some of the important properties of water as a solvent 
are reviewed in Box 4.1. 
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